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Abstract: We investigate logarithmically enhanced electromagnetic corrections of all an- 
gnlar observables in inclusive B — )• We present analytical results, which are 

supplemented by a dedicated Monte Carlo study on the treatment of collinear photons in 
order to determine the size of the electromagnetic logarithms. We then give the Standard 
Model predictions of all observables, considering all available NNLO QCD, NLO QED and 
power corrections, and investigate their sensitivity to New Physics. Since the structure 
of the double differential decay rate is modified in the presence of QED corrections, we 
also propose new observables which vanish if only QCD corrections are taken into account. 
Moreover, we study the experimental sensitivity to these new observables at Belle II. 
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1 Introduction 

By now the LHC experiment has not discovered any new degrees of freedom beyond the 
Standard Model (SM). In particular, the measurements of the LHCb experiment and the 
B-physics experiments of ATLAS and CMS have confirmed the simple Cabibbo-Kobayashi- 
Maskawa (CKM) theory of the SM [1-3]. This corresponds to the general result of the B- 
factories [4, 5] and of the Tevatron B-physics experiments [6] which have not indicated any 
sizable discrepancy from SM predictions in the B-meson sector (for reviews see refs. [7-9]). 
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However, recently the first measurement of new angular observables in the exclusive 
decay B —?■ K* ii~ has shown a kind of anomaly [10]. Due to the large hadronic uncer¬ 
tainties it is not clear if this anomaly is a first sign for new physics beyond the SM, or a 
consequence of hadronic power corrections; but of course, it could turn out to just be a sta¬ 
tistical fluctuation (see e.g. refs. [11-19]. The LHCb analysis based on the 3 fb“^ dataset is 
eagerly awaited to clarify the situation. More recently, another slight discrepancy occurred. 
The ratio Rk = BR(H+ —)• K~^fj ,~)—)• K~^e'^e~) in the low-g^ region being 
the di-lepton invariant mass) has been measured by LHCb showing a 2.6(T deviation from 
the SM prediction [20]. In contrast to the anomaly in the rare decay B —)• K* which 

is affected by unknown power corrections, the ratio Rk is theoretically rather clean. This 
might be a sign for lepton non-universality (see e.g. refs. [21-30]). 

The inclusive decay mode B —)■ is one of the most important, theoretically 

clean modes of the indirect search for new physics via flavour observables (for a review 
and updates see refs. [31-33]); especially it allows for a nontrivial crosscheck of the recent 
LHCb data on the exclusive mode [17, 34]. 

The observables within this inclusive mode are dominated by perturbative contribu¬ 
tions if the cc resonances that show up as large peaks in the dilepton invariant mass 
spectrum are removed by appropriate kinematic cuts - leading to so-called ‘perturbative 
di-lepton invariant mass windows’, namely the low di-lepton mass region 1 GeV^ < s = 
< 6GeV^, and also the high dilepton mass region with > 14.4 GeV^ (or 
> 14.2GeV^). In these regions a theoretical precision of order 10% is in principle 
possible. 

By now the branching fraction has been measured by Belle and BaBar using the sum- 
of-exclusive technique only. The latest published measurement of Belle [35] is based on a 
sample of 152 x 10® BB events only, which corresponds to less than 30% of the dataset 
available at the end of the Belle experiment. Babar has just recently presented an analysis 
based on the whole dataset of Babar using a sample of 471 x 10® BB events [36] which 
updated the former analysis of 2004 [37]. 

In the low- and high-dilepton invariant mass region the weighted averages of the ex¬ 
perimental results read 

B{B )“P = (1.58 ± 0.37) x 10“® , (1.1) 

B{B = (0.48 ± 0.10) X 10"® . (1.2) 

All the measurements are still dominated by the statistical error. The expectation is that 
the final word of the present B factories leads to an experimental accuracy of 15 — 20%. 

In addition, Belle has presented a first measurement of the forward-backward asym¬ 
metry [38] and Babar a measurement of the CP violation in this channel [36] . 

The super flavour factory Belle H at KEK will accumulate two orders of magnitude 
larger data samples [39]. Such data will push experimental precision to its limit. This is the 
main motivation for the present study to decrease the theoretical uncertainties accordingly. 

The theoretical precision has already reached a highly sophisticated level. Let us briefly 
review the previous analyses. 
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• Within the inclusive decay mode B —?• the dominating perturbative QCD 

contributions are calculated up to NNLL precision. The complete NLL QCD contri¬ 
butions have been presented [40, 41]. For the NNLL calculation, many components 
were taken over from the NLL calculation of the B —?■ mode. The additional 

components for the NNLL QCD precision have been calculated in refs. [42-54]. 


• If only the leading operator of the electroweak hamiltonian is considered, one is led to 
a local operator product expansion (OPE). In this case, the leading hadronic power 
corrections in the decay B —)• Xs£~^i~ scale with l/m^ and l/mf only and have 
already been analysed [55-60]. Power correction that scale with 1/m^ [61] have also 
been considered. They can be calculated quite analogously to those in the decay 
B —7- Xg'y. A systematic and careful analysis of hadronic power corrections including 
all relevant operators has been performed in the case of the decay B —?■ A^y [62]. 
Such analysis goes beyond the local OPE. An additional uncertainty of ±5% has been 
identified. The analysis in the case of —)• Xsi'^i~ is fully analogous and work in 

progress. There is no reason to expect any large deviation from the B —)• Xg'y result. 

In the high-g^ region, one encounters the breakdown of the heavy-mass expansion 
(HME) at the end point of the dilepton mass spectrum: Whereas the partonic contri¬ 
bution vanishes, the 1/m^ and 1/m^ corrections tend towards a finite, non-zero value. 
Contrary to the end-point region of the photon-energy spectrum in the B —> A^y 
decay, no partial all-order resummation into a shape-function is possible. However, 
for the integrated high-g^ spectrum an effective expansion is found in inverse pow¬ 
ers of = mb X (1 — ^Smin) instead of mb [63, 64]. The expansion converges 
less rapidly, and the convergence behaviour depends on the lower dilepton-mass cut 
Smin = [52]. 

The large theoretical uncertainties could be significantly reduced by normalizing the 
B —7- Xsi~^i~ decay rate to the semileptonic B —)• X^in decay rate with the same s 
cut [60]: 



dr{B Xsi^£-) 
ds 



dr(H° ^ xji^) 

ds 


(1.3) 


Eor example, the uncertainty due to the dominating 1/m^ term could be reduced 
from 19% to 9% [65]. 


• In the inclusive decay B —?■ Xsi~^i~, the hadronic and dilepton invariant masses are 
independent kinematical quantities. A hadronic invariant-mass cut is imposed in the 
experiments. The high-dilepton-mass region is not affected by this cut, but in the 
low-dilepton mass region the kinematics with a jet-like A^ and m\ < m^A implies 
the relevance of the shape function. A recent analysis in soft-collinear effective theory 
(SCET) shows that by using the universality of the shape function, a 10 — 30% reduc¬ 
tion in the dilepton-mass spectrum can be accurately computed. Nevertheless, the 
effects of subleading shape functions lead to an additional uncertainty of 5% [66, 67]. 
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A more recent analysis [68] estimates the uncertainties due to subleading shape func¬ 
tions more conservatively. By scanning over a range of models of these functions, 
one finds corrections in the rates relative to the leading-order result to be between 
— 10% to -|-10% with equally large uncertainties. In the future it may be possible to 
decrease such uncertainties significantly by constraining both the leading and sub¬ 
leading shape functions using the combined B —)• Xg'j, B —?• Xuiu and B —?■ Xsi~^i~ 
data [68, 69]. In [70], B Xsi+r in the presence of a cut on mxs was analysed by 
performing the matching from QCD onto SCET at NNLO, and a prediction of the 
zero of the forward-backward asymmetry in this semi-inclusive channel was provided. 

• As already discussed, the cc resonances can be removed by making appropriate kine¬ 
matic cuts in the invariant mass spectrum. However, nonperturbative contributions 
away from the resonances within the perturbative windows are also important. In 
the KS approach [71, 72] one absorbs factorizable long-distance charm rescattering 
effects (in which the B —)> Xgccr transition can be factorized into the product of sb 
and cc color-singlet currents) into the matrix element of the leading semileptonic op¬ 
erator Og. Following the inclusion of nonperturbative corrections scaling with 1 /, 
the KS approach avoids double-counting. For the integrated branching fractions one 
finds an increase of (1 — 2)% in the low-q'^ region due to the KS effect, whereas in 
the high-g^ region there is a decrease of ~ 10%, which is still below the uncertainty 
due to the l/rufe corrections. 

• The integrated branching fraction is dominated by this resonance background which 
exceeds the nonresonant charm-loop contribution by two orders of magnitude. This 
feature should not be misinterpreted as a striking failure of global parton-hadron 
duality [73], which postulates that the sum over the hadronic final states, including 
resonances, should be well approximated by a quark-level calculation [74]. Crucially, 
the charm-resonance contributions to the decay B —>■ Xsi^i~ are expressed in terms 
of a phase-space integral over the absolute square of a correlator. For such a quantity 
global quark-hadron duality is not expected to hold. Nevertheless, local quark-hadron 
duality (which, of course, also implies global duality) may be reestablished by resum¬ 
ming Coulomb-like interactions [73]. 

• Also electromagnetic perturbative corrections were calculated: NLL quantum elec¬ 
trodynamics (QED) two-loop corrections to the Wilson coefficients are of 0(2%) [54]. 
In the QED one-loop corrections to matrix elements, large collinear logarithms of the 
form \og{ml/m‘j) survive integration over phase space if only a restricted part of 
the dilepton mass spectrum is considered. These collinear logarithms add another 
contribution of order -|-2% in the low-g^ region of the dilepton mass spectrum in 
B —)■ XsfJ.^fJ.~ [75]. For the high-( 7 ^ region, one finds —8% [65]. 

Based on all these scientific efforts of various groups, the latest theoretical predictions have 

been presented in ref. [65]. 
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In the present manuscript, we make the effort to provide all missing relevant pertur¬ 
bative contributions to all independent observables in the decay B —)• As it is 

well-known, the angular decomposition of this inclusive decay rate provides three indepen¬ 
dent observables, Ht, Ha, Hl from which one can extract the short-distance electroweak 
Wilson coefficients that test for possible new physics [76]: 

^^2 ^ 8 “ z‘^)HL{q^) + ■ (1.4) 

Here, z = cos0, where 0 is the angle between the and B meson three momenta in the 
di-lepton rest frame, Ha is equivalent to the forward-backward asymmetry [77], and the 
spectrum is given by Ht + Ht- The observables dominantly depend on the effective 
Wilson coefficients corresponding to the operators 07 , 0 g, and Oig. 

The paper is organized as follows. In section 2 we define the observables which we 
consider in the present analysis. In section 3 the derivation of the log-enhanced terms is 
presented. Master formulae for our observables are given in section 4, our phenomenological 
results in section 5. We briefly discuss the new physics sensitivity of our observables in 
section 6. Finally we explore the precise connection between experimental and theoretical 
quantities using Monte Carlo techniques in section 7. The latter analysis updates, and 
in parts supersedes, our previous statements in ref. [78]. We conclude in section 8. In 
the appendices we collect various functions that arise in the computation of QED and 
QCD corrections to the observables (appendix A), as well as formulas that parametrise the 
observables in terms of ratios of high-scale Wilson coefficients (appendix B). 


2 Definition of the observables 


The z dependence of the double differential decay distribution presented in eq. (1.4) is exact 
to all orders in QCD because it is controlled by the square of the leptonic current. The 
inclusion of QED bremsstrahlung modifies the simple second order polynomial structure 
and replaces it with a complicated analytical z dependence (see eqs. (3.28)-(3.33)). In 
particular this implies that, as long as QED effects are observably large, a simple fit to a 
quadratic polynomial will introduce non-negligible distorsions in the comparison between 
theory and experiment. In this section we explain the procedure that we adopt to construct 
various differential distributions and suggest that experimental analyses follow the same 
prescriptions. 

The extraction of multiple differential distributions from eq. (1.4) is phenomenologi¬ 
cally important because the various observables have different functional dependence on the 
Wilson coefficients. Eor instance, at next-to-leading order in QCD and without including 
any QED effect, the three Hi defined in eq. (1.4) are given by [76]: 


Hriq^) 

HL{q^) 


487r3 

Gimi\ytsyth? 

487r^ 


2s(l — s)^ 111^9 + ^^71 + |Cio|^ 


(1 — s)^ 1^9-1-2(771 + |(7io|^ 


( 2 . 1 ) 

( 2 . 2 ) 
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HA{q^) 


(2.3) 


G\ 


ml 


\Vt*sVti,Y 


487r3 


(-4s) (1 - sf Re 


<^10(^9 + T Ct ) 


We decided to preserve the natural definitions of the differential decay width dF/ dq^ 
and of the forward-backward asymmetry dAYB/dq^'. 


dq^ 

d^FB 

dq^ 



d?T 

dq^dz 

d^r 

dq^dz 


dz , 

sign(z)d 2 : , 


(2.4) 

(2.5) 


with the understanding that ^fb does not coincide with the coefficient of the linear term 
in 2 ; in the Taylor expansion of d^T/dc^dz. 

We extract other single-differential distributions by projecting the double-differential 
rate onto various Legendre polynomials, Pn{z). These polynomials are orthogonal in the 
[—1,1] interval and are, therefore, ideally suited as angular projectors. In order to make 
connection with the existing literature we choose the first two projections in such a way 
to reproduce Ht and Hi in the limit of no QED radiation. For the higher order terms we 
simply adopt the corresponding Legendre polynomials. The observables are defined as 


Hi{q^) = (2.6) 

and the weights we use are: 

Wt=‘^ Pq{z) + y P2{z) , 

Wl = \po{z)-^-Ip2{z), 

4 

= gSign(z) . 

Note that Wt + Wl = Po{z) = 1 implying that the relation dV/dq^ = Ht + Hl is exact. 

The unnormalized (defined in eq. (2.5)) forward-backward asymmetry receives contri¬ 
butions form all odd powers of 2 ; in the Taylor expansion of the double differential rate and 
is given by 


VFs = Psiz) , 

W 4 = P 4 (z), (2.7) 


dA^B 

dq^ 


\ HA{q^) ■ 


( 2 . 8 ) 


In the literature the normalized differential and integrated forward-backward asymmetries 
are often considered (see for instance ref. [65]): 


d^FB 

dq^ 


d^T 


j d^T 


3 HA{q^) 

AHT{q^)+HL{q^) ’ 


(2.9) 
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Figure 1. Kinematics of collinear photon radiation. The collinear photon is radiated off (left 
panel) and l~ (right panel), respectively. The crossed grey circles denote an operator insertion 
from the effective weak Hamiltonian. The arrows indicate momentum rather than fermion flow, x 
denotes the momentum fraction of the collinear photon. 


^FB[qm,qM] 



fill 


3 

4' 


dq^HA{q^) 


fill 


dq^ [Hriq^) + Hiiq^)] 


( 2 . 10 ) 


The new observables H'^ and H 4 (obtained by employing the weights and W 4 ) 
vanish exactly in the limit of no QED radiation but are still potentially important for 
phenomenology because of their non trivial dependence on the Wilson coefficients. We 
find that projections with even higher Legendre polynomials are suppressed and will not 
be considered further. 

Note that the expected statistical experimental uncertainties (at a given luminosity) are 
well understood in the total width {Ht + H^) and forward-backward asymmetry (3/4 Ha) 
cases. On the other hand, Ht, Hl, H 3 and are obtained by projecting the double 
differential rate with weights that (especially for W 3 and IL 4 ) are essentially arbitrary. As 
a consequence a simple rescaling of these weights implies a corresponding rescaling of the 
central values we find. In section 6 we show how to use the squared weights {Wf) to assess 
the expected Belle II reach for each of these observables. 

The experimental procedure that we recommend is to use the weights Wj to extract 
single-differential distributions and to refrain from attempting polynomial fits to the data. 


3 Log-enhanced QED corrections to the double differential decay rate 

In this section we work out the formulas for the logarithmically enhanced electromagnetic 
corrections of the double differential decay rate d^T/{dq^ dz). The operators and Wilson 
coefficients of the effective weak Hamiltonian are the same as in [65, 75]. The kinematics 
can be inferred from figure I. 

Let us first consider the case without photon radiation. The momenta of the quarks 
are labelled ph and ps, respectively. The momenta of and i~ are denoted by pi and p 2 , 
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respectively. From momentum conservation we arrive at Pb = Pi P 2 + Ps- We define the 
invariants 

‘^PiPj 


Sij — 


m 


2 ’ 


Moreover, we define 


yi = 


b 

2 El 


i G { 1 , 2 , 5,6} . 

2 E 2 

y 2 = 


(3.1) 


(3.2) 


mb mb 

where Ei [i = 1 , 2 ) is the zero-component of pi when evaluated in the rest-frame of the 
decaying 6 -quark. From momentum conservation and by treating all hnal-state particles 
as massless, we obtain the relation yi + y 2 = ^ + •512 • This relation also implies 


•Sis — 1 — 2/2 , 

S2s = ^ — yi 1 


•Sift = 2/1) 
S2b = 2/2 • 


■Ssfe = 1 — 'S12 , 


(3.3) 


For the double differential decay rate we also need 2 ; = cos 9, where 9 is the angle between 
the 6 -quark and the positively charged lepton in the centre-of-mass system (c.m.s.) of the 
final-state lepton pair. Hence 

z = cos0=^-i, (3.4) 

\Pi\ \Pb\ 

where all primed momenta are taken in the c.m.s. of the final-state lepton pair. It turns 
out that z is simply given by [56] 

2 = . (3.5) 


1 - Sl2 

At this point we stress that the LHS of this equation is evaluated in the lepton c.m.s., 
whereas its RHS is evaluated in the rest-frame of the decaying 6 -quark. The connection 
between the angle 9 and the leponic energy asymmetry has already been emphasized in [56] . 

We now switch on QED and consider the radiation of a collinear photon off a lepton 
leg as shown in figure 1. The momentum of the positively (negatively) charged lepton is 
denoted by pi (^ 2 ) before radiation and by (p_) thereafter. If the positively charged 
lepton radiates the photon (left panel of hgure 1 ), its momentum after radiation is 
given by = xpi, where x denotes the momentum fraction of the collinear photon and 
X = 1 — x. In this case, the momentum of the negatively charged lepton remains unchanged 
and hence we have p- = p 2 - If the negatively charged lepton radiates the photon (right 
panel of figure I), we obviously have p- = xp 2 and = pi. In analogy to eq. (3.2), we 
define 


y± = 


2 E± 

mb 


(3.6) 


where E± is the zero-component of p ±, again evaluated in the rest-frame of the decaying 
6 -quark. We will also need the definition 


S+_ = X Sl2 


(3.7) 


As already discussed in refs. [65, 75], the logarithmically enhanced contributions stem¬ 
ming from collinear photon radiation are evaluated by 

^ bcoii _ d Fqqji 2 d Fqqii 3 


ds dz 


ds dz 


ds dz 


(3.8) 
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where we stay differential in the double invariant s_|_= (p+ +p_)^ = x si 2 and the triple 

invariant S 12 = (p+ + P- +^ 7 )^ = (pi +P 2 Y 1 respectively. We first consider the case of 
the triple invariant, where the formulae look exactly the same as in the case without QED, 
since we can lump the lepton and the collinear photon. We therefore arrive at 


drcoii,3 = PFdsi2dyidy2dx5{l +si2-yi-y2) (si2, yi, ^ 2 ) 

x6»(yi) 0(1 - yi) 0 ( 2 / 2 ) 0(1 - 2 / 2 ) 0 (si 2 ) 0(1 - 512 ) 0 {x) 0(1 - x). 


(3.9) 


Here 0 denotes the heaviside step-function, PF is the pre-factor 


PF = 


Glmb\VtbVt*s? 

327r3 


(3.10) 


and f\^\x) is the mass-regularised splitting function for collinear photon radiation of 
which we only keep the logarithmically enhanced part (5e = Q;e/(4vr)), 



4ae- 


[1 + (1 - x)2] 


In 



(3.11) 


The squared matrix elements \A\‘^ for the different operators read 
(si2,2/1, 2/2) = - - [(1 - 2/2) 2/1 + (1 - 2/1) 2/2] , 

•512 

l^lra (512,2/1,2/2) = ‘^rnf (1 - S12) , 

l^lgg (512,2/1,2/2) = 4 mf{l - 2/i)(l - 2/2) + 2 mfe S12 (1 - S12), 
l-^lno (512,2/1,2/2) = ^rnf (2/1 - 2/2), 
l-^lgio (512,2/1,2/2) = 2 mf S12 (2/1 - 2/2), 

1^1 29(512,2/1,2/2) 0^6/2(512) l-^lgg (512,2/1, 2/2) , 

1^1 27(512,2/1,2/2) — Cle/ 2 ( 5 l 2 ) |- 4 .| 7 g (Sl 2 , 2 / 1 , 2/2) , 

\A\I2 (512,2/1,2/2) = 5e 1/2(512)!^ I^Igg (512,2/1,2/2) , 

1^1210 (512,2/1,2/2) = 5 e f 2 {si 2 ) l^lgio ( 512 , 2 / 1 , 2 / 2 ) • ( 3 . 12 ) 


The function / 2 (si 2 ) denotes the one-loop matrix element of P 2 and is given by 

In — m, when 2/c < 1 

l-yfl-Vc ’ 

2 arctan ^ , when 2/c > 1 




Vi/c-i ’ 


(3.13) 


with yc = 4mc/(m^si2). 72 ( 512 ) is a complex function and therefore the \ A\ 2 j with j / 2 
are complex. However, after taking into account the Wilson coefficients and adding the 
appropriate complex conjugate expression, the double differential rate turns out to be real, 
see eq. (3.28). 

Let us now come back to the evaluation of (3.9). After integrating over the J-function 
and changing variables according to eq. (3.5) we arrive at 
1 

= 2 PF [dx6^\x) 
dsi 2 dz J ^ 


( 512 , 


1 + 512 1 — 512 1 + 512 1 — 512 

- 2 ;,-- 2 : 

2 2 2 2 


0 
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Pi +P2 




Figure 2. Boosts that have to be performed in order to translate the b rest-frame into the leptonic 
c.m.s. Left panel: without QED. Middle panel: Radiation off i~. Right panel: Radiation off 


X 0(1 - z) 9(1 + z) 0 ( 512 ) 0(1 - S 12 ) ■ (3.14) 

The factor of two stems from the fact that both diagrams in figure 1 are relevant. Note 
that the integral in eq. (3.14) is divergent at x = 0. However, eq. (3.8) is well-behaved 
once all expressions on its RHS are plugged in. 

We now turn our attention to the more complicated case of the double invariant 
coll, 2 /{ds dz), and hrst address radiation off £~. As can be seen from the middle panel 
of figure 2, the boost A from the 6-quark rest-frame into the leptonic c.m.s. is determined 
by 

PI + XP2^0. (3.15) 


After the boost, we compute z = cos 9 via 

^ ^ p+ • Pfe ^ xy2 - yi 

l^+l \Pb\ ^ (y^ + xy2f - 4X512 


(3.16) 


Again, the primed momenta are evaluated in the lepton c.m.s., whereas the RHS of the 
equation is evaluated in the rest-frame of the 6-quark. The differential decay width reads 

^^rcoii,2 = PFdsi2dyidy2dx5{l +si2-yi-y2) (si2, yi, ^ 2 ) 

x0(yi) 0(1 - yi) 0 ( 2 / 2 ) 0(1 - y 2 ) 9 {si2 ) 0(1 - S 12 ) 0(x) 0(1 - x). (3.17) 


We hrst eliminate yi by integrating over the (5-function. Subsequently, we eliminate 2/2 in 
favour of z according to eq. (3.16). This transformation reads 

(±), N (1 + S12)(2 - X - XZ^) ± 2z\/l - xJ (1 - 512)^X - 512x2(1 - Z^) 

= - x^(l-.^)+4x - 

It turns out that this in an injective mapping only for 512 < x. For 512 > x we have 
to subdivide the 2 / 2 -interval into two pieces, so that we get a total of three contributions. 
After the additional variable substitution 512 = s-\ _/x they read 


(3.18) 


“ coll,2; 1 

ds-i_ dz 


PF 



0 


\ d A+) 

X 5z^2 



(■S 12 , 1 -b 512 


02 , 2 / 2 ) 


02 = 0^+^^) 

S 12 = s+-/x 
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X0(1 — z) 0(1 + z) 0(s_| _) 0(1 — s_|_) , 


(3.19) 


^ 2 -p£ 

“ coll,2; 2/3 

ds-i_dz 


where 


= ±PF 


dx 


ft\x) 


9/2 


1^1 (si 2 , 1 + S 12 - 2 / 2 , 2 / 2 ) 


xO(-z) 0(1 + z) 0(s+_) 0(1 - s+_) , 

1 - s+ 


(=t) / \ 

02 = 02 (^) 
S 12 = s+_/x 

(3.20) 


- /n '2^ —" • 

1 T \/(l - ^2)s^_ 

Once the photon is radiated off £~^, we apply very similar steps. As can be seen from 
the right panel of figure 2, the boost A, from the 6-quark rest-frame into the leptonic c.m.s., 
is determined by 


xpi + P2 

After the boost, we compute z = cos 9 by 


z = 


P+ • Pb 
\P+\ \Pb\ 


0 . 


02 - xyi 


xyi + 2 / 2 )^ - 4xsi2 


(3.22) 


(3.23) 


We now eliminate 2/2 by integrating over the d-function. Subsequently, we eliminate 2/1 in 
favour of z according to eq. (3.23). This transformation reads 


(±), N (1 + ■S12)(2 - X - XZ^) ± 2z\/l - xJ (1 - S12)2x - Sl 2 x‘^{l - Z^) 
01 (z) = - 


x2(l-z2)+4a: ' 

As mentioned before, this is an injective mapping only for S 12 < x. For S 12 > x we have 
to subdivide the 2 / 1 -interval into two pieces, so that in this case we also get a total of three 
contributions. After the additional variable substitution S 12 = s_|_ /x they read 

‘Pp.iua ‘^7^ fP(^) [S , 


1^1 (Sl2,01,l + S 12 - 01 ) 


X0(1 — z) 0(1 + z) 0(s-|_) 0(1 — S-I_) , 


01 = 0 ^”/) 
S 12 = s+-/x 
(3.25) 


j2-p£+ 

“ ^ coll,2; 2/3 

ds_|_ dz 


= 


dx 






|-21P(S12,01,1 + S12 -0l) 


1 -W+- 

X0(z) 0(1 — z) 0(s_|_) 0(1 — s_|_) . 

The total contribution in case of the double invariant is now obtained by 

d^rcoii,2 

i=l 


ds^ _ dz 


E 


V/ 2 p^+ j 2 p£ 

^ ^ coll,2; i coll,2; i 

ds_i_ dz ds_|_ dz 


01 = 0 /( 2 ;) 

Sl2 = S+-/X 

(3.26) 


(3.27) 
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We finally identify S 12 = s in eq. (3.14) and s_|_= s in eq. (3.27) and plug everything 

into eq. (3.8). This leads us to the following expression for the logarithmically enhanced 
collinear decay width 


d^rcoii 
ds dz 


G ~ 1 (^b\ \ \r^ \2 ^ i 1/^ 12 

-ln( 32 ) <( ICgl '(s,2;) + |Cio| ^99 (■SW) 


Km}. 


Tag ? 77 ™^('S, z) + Qfe R-e 

+Re [C9Ci*o] ^) + 52 Re [(C 2 + CpCi) Cf* z) 


/^eff 

O 7 Og 


+ Oe Re 


/^eff 

O 7 0;LQ 


f{‘ioT.=) 


“ho^e 


{C2 + CFCi)CUi"\s.z) 


+ (Ue Re 


\C2 + CFCi)Cl^i^^\s,z) 


+al{C2 + CFC^_f itr\s,z) 


(3.28) 


where we assumed that the Wilson coefficients Ci and C 2 are real, and we neglected 
contributions from the penguin operators P 3-6 due to their small Wilson coefficients. The 
functions ^^-^“^( 5 , 2 ;) are given by 


,(=„), , 64p,(s,j)^l„(y3-y: 

S77 I'S) Z) — 


- 1 


+ 


+ 


{z‘^ — 1)^ Vs + z"^ — 1 
64pVs,z) 


+ 




3T-1 


+ 


6 Azp 2 {s,z) 

s{z‘^ — 1)3 

16 p 4 (s, z) ln(s) 


s{z^ - 1)3 (s (z2 _ 1) + 1)3/2 s(^2 _ 1)3 

ApVs,z) 16 (s-l) 2 p 6 (s,z) ln(^M^ 




3s {z^ - if {s (z^ - 1) + 1) 

16 szp 7 is,z) In(^) Asp 8 {s,z) ln(s) 


+ 


+ 


Pio{s,z) 


^iT\s,z) = - 


(z2 - if (z2 - 1)4 

8s3/2p9(s,z) ln(^/^-^^- 

(z2 — 1)4 (s _|_ 2:2 — 1)5/2 3|-^2 _ 1)3 ^2 _ 1)2 

+4(s-l)q«^ + ,-^= + l)ln/?|^) , 

64zpn(s,4 In(iTl) 32pi2(s,z) lii(s) 


(z^ - 1)3 (z2 _ 1)3 

8 pi 3 (s,z) , MV~spiVs,z)\n(^,J^ - 


- 1 


(Z2-1)"(S + Z2-1) 


+ 


(z2 - 1)3 (s + z2 - 1) 


3/2 


(».*-) = - 


64pi5(s,2) sigti(4)ln( WVi ** —) 




(z2 - 1)2 y/s (z2 - 1) + 1 


(3.29) 


(3.30) 
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t(em) 

?910 


(s,z) 


64pi6(s, z) sign(z) 


32 


+ 


(^2 - 1)3 (s + z2 - 1)3/2 

zpi7{s,z) ln(i(Vs + l)Vl-^2^ 8(V^-l)"zpi8(s,2) 


(z2 - 1)3 

64 s z (9 sz2 + 7s + 42;2 — 4) ln(s) 


+ 


(2;2 — ly (s -I- 2;2 _ 1) 
— 32(s — 1)^ z ln(l — -v/s) , 


(3.32) 


+ 


(z2 - 1)3 

32 s zPi 9 (s, z) ln(s) 16 szp 2 ois, z) ln(i (+ 1) Vl - 
(z2 - 1)^ ^ (z2 _ 1)4 


(z2-l)3 (s + ^2_l)2 
— 16(s — 1)^ s z ln(l —-v/s) . 


(z2 — 1)^ (s + ^2 — 1) 


5/2 


(3.33) 


The pi{s,z) are polynomials in s and z and are given in appendix A. In case of negative 
or complex arguments, the logarithms and square-roots are defined as 

^ = ^e'/ 2 -gW , 

ln( 2 ;) = ln|z| + i arg(z) , 

arg( 2 ;) G (—vr, vr] . (3.34) 


The functions z) with X = 2, 7, 9,10 cannot be computed analytically since the 

squared matrix elements (see eq. (3.12)) are complicated functions of si 2 . We therefore 
refrain from presenting their explicit expressions. They can easily be computed numerically 
by applying the steps outlined above. 

A strong cross-check is done if we weight the ^^j'^\s,z) by unity or by sign(z) and 
subsequently integrating over 2 . After proper normalisation one obtains the functions 
from [65, 75]. Note that this cross-check is non-trivial due to the fact that in 
our former work we computed the w|j™^(s) in a different way: Since there was no need 
to introduce the variable z we performed the calculation entirely in terms of the rescaled 
energies Moreover, there was more freedom in choosing the order of integrations since we 
were not forced to perform the x-integration immediately after that over the 5-function. 
These two simplifications led to significantly simpler variable substitutions and shorter 
expressions. With the ability to reproduce them by the more complicated calculation can 
therefore be regarded as a non-trivial cross-check. 

4 Master formulas for the observables 

We start again from the double differential decay width 

y 2 p o 

'jyyyp ^ g [(^ + + 2(i - z‘^)HL{q^)], (4.1) 
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where z = cos 9 and 9 is the angle between the and the B meson three momenta in the 
di-lepton rest frame. This formula is modified once QED corrections are taken into account 
(see sections 2 and 3) due to the appearance of higher powers of 2 . As stated in section 2, 
we project out the Hj (/ = T, A, L) by eq. (2.6) in this case. Then the Hj are functions 
of the dilepton-invariant mass but obviously not of z. Ha is proportional to the 

lepton forward-backward asymmetry; the g^-spectrum is given by Hx -|- , 




-1 


(fV 

dz dq^ 

a 


= Hxiq^) + HL{q^) , 


dApB ^ d‘^r . 3 2 ^ 


Each of the Hj can be expressed as follows {s = 0^/"ml poie)- 




(4.2) 

(4.3) 


(4.4) 


where the dimensionless functions ‘^^^(s) include both perturbative and non-perturbative 
contributions. Moreover, we normalise the observables to the inclusive semi-leptonic b —)• 
AcCP decay. However, the normalisation proceeds in such a way that we insert the perturba¬ 
tive expansion of the inclusive semi-leptonic b —)• X^eu decay (including power-corrections), 
and also use the ratio [54, 79] 


Vub ^ T{B^ X,ev) 
Vcb T{B ^ Xucu) ’ 


(4.5) 


which was recently reanalysed in [80]. We therefore use C = 0.574±0.019 (see also table 1). 
Consequently, our expression of the normalised angular observables Hj reads 


Hi 


B{B —)• Acei^)exp 


VtlVtb " 4 4>|,(g) 

Vcb C 


(4.6) 


where <!>„ is defined by [75] 

r(H —)• X^eu) 


G'^pm 


5 

b,pole 


1927r3 




(4.7) 


The expansion of is given by 


$„ = 1 -h -F K 


12 / _ix 


+ 


¥;( 2 )+ 2 /?fV(i)ln(^ 

\mb 




Ai 9 A 2 


+ 0{a^,K^,asK,asA^/ml,A^/ml) , 
50 Svr^ 

^ ~ Y~ ^ ' 


(fV) = rih 


2048C3 16987 3407r2 

+ -I + ni 


9 


54 


81 


[256(3 1009 3087r2 


V 9 


27 


81 


41848C3 5787r^ 1044807r2 1571095 848 2 , , x 


2 m^ 


2 m^ 


(4.8) 
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As explained in detail in [75], a consistent perturbative expansion in inclusive B —)• 
in the presence of QED corrections is done in 5s = Q!s(/ib)/(47r) and k = ae{nb)/cis{fJ'b)- 
We will also briefly sketch the structure of this expansion later below. 

In the above equation, the 0{k) is taken from [81], with r] = as{fJ-o)/as{^ib)- There 
also exist QED corrections at 0(5 sk) which could be computed in principle. However, 
they are not logarithmically enhanced since the fully integrated B —?■ rate is an 

infrared safe observable with respect to collinear photon radiation. We therefore neglect 
this contribution, but include it lateron in the quantity TZ{so), where QED logs will be 
present in the normalisation. 

The two-loop correction of 0(5^) was taken from [82]. Here, Uh and ni are the numbers 
of heavy and light quark flavours, respectively, and /3q = 23/3 is the one-loop QCD 

/^-function for five active flavours. The coefficients Ai and A 2 in the power-suppressed 
terms represent the matrix element of the kinetic energy and magnetic moment operator, 
respectively, and are defined as 


Ai = {B\h{iDfh\B)/{2MB), 

A 2 = -{B\hiaf^^G^,h\B)/{12MB) ~ - M|) . (4.9) 

As far as the quantity ^£i{s) is concerned, we expand it in the following way in terms of 
products of the low-scale Wilson coefficients and various functions arising from the matrix 
elements. 


i<j 


Cfif,b)Cf*{fib)H!^if^b,s) 


(4.10) 


where Cf^{^b) / Ci{^b) only for i = 7, 8. Here i and j run over all operators of eqs. (15) 
and (16) in [75]. Their low-scale Wilson coefficients are also given explicitly (analytically 
and numerically) in that paper. Eor I = T,L the functions s) are given by 




^ = < 


for i = j ^ 


iV=7,9,10 


E 2M^M^* + (M/Mf + MfM]*] Sjg + ARI , Iot i ^ j . 


(4.11) 




iV=7,9,10 

\ 

For I = A the formula is simpler, 




0 , for i = j , 

E + M/OMf + AHf. , ioi i ^ 3 . (4-12) 


N=7,9 


The coefficients Mf are listed in table 6 of [75]. The building blocks have the following 
structure. 


^NM — 


^ nm { s ) |l + 8as w^M^/('S) + 16ag w^]v4-^/(s)| 



(4.13) 
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From (4.11) and (4.12) we see that the possible combinations of indices are NM = 
77, 79, 99, 1010 for I = T,L and NM = 710, 910 for I = A. Moreover, we have Sgg = S'foio 
for I = T,L. Explicitly, the phase-space factors read 


a^^(s) = 8(1 - s)Vs, 

= 4(l-sf, 

= -8(1 - sf , 

o-79(s) = 8(1 - s)^ , 

^79{s) = 4(l-s)2. 

= -4s(l-s)^ 

o-^9(s) = 2 s(l-s)^ 

<^99(«) = (1 - • 

(4.14) 

The one-loop QCD functions 

extracted from [49] and have already 


f<y\ 

been given in [76]. The two-loop functions j{s) have so far only been available for the 

q^-spectrum [83-86], but not for the double differential rate. Due to a recent calculation 

of the double differential rate of the inclusive semi-leptonic b —>• decay at two loops 

in QCD [87], they can be extracted for NM = 99, 1010 and 1 = T,L as well as for 

NM = 910 and I = A. The data to extract these functions was kindly provided by 

the authors of [87, 88] and we can therefore present them here for the first time. All 
f 1 

functions j{s) are rather lengthy and we therefore relegate their explicit expressions 
to appendix A. 

The functions xInm^^) ~ accompany the non-perturbative 0{ls!^Qjg/m^) 

corrections can be obtained from [56] (see also [55, 58]) and were previously computed 
in [76]. We confirm their expressions, 



1 = (1 - s)(5s-F3), 

3s 

Xi,77is) = ^(s-l)(3s + 13), 


Xg79(®( 

I = 4(1-s)2, 

Xl',79('S) = 2(1 — s)^ , 


xf,99(®; 

' = “ s)(3s -k5), 

Xi', 99('5) = g(l — s)(13s -1- 3), 


Xl^71o(®; 

1 = — — (3s^ -|- 2s 4- 3) , 

0 

Xg9io('5) — “3® (8'^^ -1- 2s -1- 3) , 

(4.15) 

X2,77{s) 

= i (3s2 -4 2s - 9) , 
s ^ 

xi,77(s) = 2 (15 s2 - 6 s - 13) , 


X2,79i^) 

= 4 (95^ - 6s - 7) , 

xf,79(s) = 2 (3s2 - 6s - 1) , 


X2,99i^) 

= s (15 s 2 - 14s - 5) , 

X2,99(®) = 2 ( —+ 10'S + 3) , 


X2,71oi^) 

= -4 (9s 2 - 10s - 7) , 

X^9io(s) = -2s (I5s^ - 14s - 9) . 

(4.16) 

The quantities can be further decomposed into 



ah!^ 


(4.17) 


Here the contributions represent finite bremsstrahlung corrections that appear at NNLO. 
They are known for the g^-spectrum (i.e. I = T + L) [47] and the forward-backward 
asymmetry (equivalent to / = A) [50], but not for the double differential rate. Hence we 
only include them for these two cases, but not for Ht and Hi separately. This is still an 
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excellent approximation since the effect of finite bremsstrahlung corrections is very small 
anyway. The explicit formulas can be found in [47, 50] and will therefore not be repeated. 

The coefficients clj comprise non-perturbative contributions and were 

calculated in ref. [61] for I = T + L and I = A. Moreover, the coefficients of the double 
differential rate can be inferred from that paper. One obtains 


T 

~ 8A2 , 


- 2 ( 

ymj 

T 

1 T 

Oj 

I <M 

lICD 

1 

II 

T 

~ 8A2 , 

C22 

- 2 ( 

T 

~ 4A2 

Cll 

“ +“^'^ 27 m 2 

T 

~ 8A2 . 

C12 

- “■'‘Smd 


for j / 1,2 , 

1 — s)^(l + 3s) F{r 

(1 — s)^(l + 3s) F{r) 


- m'^* + - M^* 

s ^ 2 ^ 


- Ml* + - Ml* 

s 2 

- MJ* + - M?* 
s ^ 2 ^ 


1 1 


for j / 1,2 


F* (r) (-Ml + -Ml\-- F{r) ( - Ml* + - Ml 


9 * 


for j / 1,2 , 


F{r)(^Ml* + ^Ml*'^ , 

~ '5)^(1 + 3s) F{r) , 

Cno = “g <^^10 ) (4-18) 

where r = g^/(4m^). The function F(r) can be found in the appendix of [75]. Moreover, 
we also include factorisable non-perturbative charm contributions which we implement by 
means of the Kriiger-Sehgal approach [71, 72]. We elaborated extensively on this approach 
and also the formulas by means of which these corrections are taken into account in ref. [65] . 
Given their length we do not repeat these formulas here but refer the inclined reader to 
refs. [65, 71, 72] for all necessary details. 

Finally, the coefficients collect the ln(m^/m^)-enhanced electromagnetic corrections 
which we calculated in section 3 for the double differential rate. Their contribution to the 
Hj can be derived from (3.28) by applying the projections given in section 2. One finds 

677 = SttsK cj 77 (s)a;^ 7 /(s) , 629 = Sa^K agg{s) ’ 

I O ~2 2 / /-N I o ~3 3 / 

679 = 8 a^K cj79(s)a;^9 /(s) , 632 = 8a,K cj 99 (s) ^2,7 («) > 

eig = 8 asK 0 -^ 9 (s) ^ 9 ^ (s) , = f 42 > 


4 j = - 4 F{r) 


M]* + -M'^ 


9 * 


4j 


1 

6 


for j / 1,2 


8A2 


C 22 = (1 - s) (3 - s) F{r) 


Cii — -bCTsK 


4 A 2 

27m2 


Ml* + ^ Ml* 


{I - s f{3 - s) F{r) 


Ml* + - Ml 
1 


1 


8A2 


C 12 = -asK^(l-s) (3-s) 


9 * 


F*{r){Ml + ^MA-^ 
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for j = 7,9, 


(4.19) 


^1010 — 699 ) 


®12 — 3 ^22 


®710 


r3„3 r 


o-no( 


,(em)/gx 

‘^27,7 Wj ) 

I 4 / 

= 3%’ 

one gets 


A 

^210 ~ 

,(em) 

A 


^110 — 


. ~2 ,^2 
I CXgK (J, 


4„A 


910 


(s) U, 


(em) /-') 
210,ylv*'' ’ 


(4.20) 


The functions a;|j™^(.s) have again been moved to appendix A. 

We consider the observables Hj (or equivalently l-Li) in the \ow-q^ region only, because 
their sensitivity to New Physics is highest in this region [76]. Besides, there are two more 
observables which we compute in the low-s region. First, there is the zero crossing 
of the forward-backward asymmetry, which we extract numerically from T-L^ by means of 
the formulas given above. Moreover, there is the branching ratio. In principle, it can be 
obtained by taking the sum of T-Lt and T-Ll- Its master formula has already been given 
in [75]. We therefore only highlight two small pieces which are available for the branching 
ratio only, but not for T-Lt and T-Ll individually. These are only the finite bremsstrahhmg 
contributions from [47] and the non-log enhanced terms of (see eq. (94) of ref. [75]). 

In the high-g^ region we consider two observables. The first one is the branching ratio, 
where we include the same terms as in the low-g^ region. As far as QED corrections are 
concerned, the functions a;gg™^(s), a;jQ™g(s), a; 77 ™^(s), and a; 7 g™^(s) (see eqs. (94) and (100) 


(102) of [75]) are valid in the entire g^-region, while the functions W 2 ^™^(s) are again 
obtained from a numerical fit. To take into account our most recent input parameters (see 
table 1), we re-did the fits and collected the results in appendix A. In addition, the two- 
loop QCD matrix element functions FJ 2 {s) and Ff 2 {s), which were originally computed 
in [52], were given explicitly only in [53]. We implement these formulas in our numerical 
code. Moreover, non-perturbative 1/m^ corrections become sizable in the high-s region. 
They were originally computed in [59] and we implement the formulas of refs. [59, 60]. The 
second observable is the ratio TZ{sq) which we have already mentioned in the introduction. 
It was proposed in [60]^ and is obtained by normalizing the B -A Xs£^i~ decay rate to 
the semileptonic -A- rate with the same cut in q^. In this way, large theoretical 

uncertainties that stem from poorly known parameters in the 1/m^ and 1/m^ power- 
corrections can be significantly reduced, as we will see in our numerical analysis in section 5. 
In terms of our perturbative quantities, it reads 

dT{B^Xsl+i-) 

.\V,lVw 


ds 


7^(so) = 


'^0 _ 

,^dT(B^ 

/ ds- 

JSQ 


XJu) 


= 4 


Vub 


ds ^„(s) 


(4.21) 


ds 


The quantity ^u{s) is known from the branching ratio. The differential <hn(s) is given by 


dr{B^^XJn) Gl\Vub\‘^m\ 


bjpole 


ds 


1927r3 


$«(s) . 


(4.22) 


^Note that we use a different pre-factor here. 
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We elaborated extensively in ref. [65] about how to obtain the 0(1,5s, Sg, 1/m^, 1/m^) 
contributions to and will therefore not repeat these formulas. We would rather 

like to describe the 0(5 sk) contribution to which we include in the present work 

and which was absent in [65]. Once the integration over s is restricted to the high-g^ 
region, the corrections of 0(5 sk) to d>u(s) contain residual terms logarithmically enhanced 
by These must be proportional to ujgg^\s}. We take into account that we 

only have one charged lepton in the final state, and that the leptonic current is y — A. 
Moreover, we average over e and fi, and arrive at 




CtsK, 


8 UgH (1 


(1 + 2 .§) ujg^^\s) 



(4.23) 


As in our previous analysis [65] we do not include electromagnetic corrections of order 0{k) 
to 4>u(s) because they are unknown. 

Let us conclude this section by a few remarks on the renormalisation schemes for the 
quark masses, as well as on the expansion in 5s and k. The pole masses of the b and c quark 
that are present in the definition of s and in several loop functions suffer from renormalon 
ambiguities [89, 90]. We therefore convert them analytically to short-distance schemes (IS 
and MS, respectively) before any numerical evaluation of the observables is carried out. In 
our numerical analysis we use the conversion formulas up to order 0(5g) [91]. As far as 
the mass of the top quark is concerned we take the pole mass as input and convert it to 
the MS scheme at order 0{af) using RunDec [92]. We also take into account electroweak 
corrections presented in eq. (31) of ref. [93], consistently to the other contributions. Turning 
our attention to the expansion in 5s and k, we observe that the amplitude has the structure 

A = K [AlO + 0!s AnLO + -^NNLO + ^(“s)] 

++ “s+ C’(as)] +0{k^), (4.24) 


and that the ratio ^££{s) /<I)„ in (4.6) has a similar structure to that of the squared amplitude 
(up to bremsstrahlung and non-perturbative corrections), 

A‘^ = A\o + ots 2 AloAnlo + o?s {A%io + ‘^‘AloAnnlo) 


+ K-^ 


TOg 2 {AnloAnnlo + .••) + 0 (ag) 


+ “s ‘^i-^NLO-^LO + -^LO-^^NLo) 

+“s + -^NNLO-^To + -^LO-^WnLo) 

+“s 2 (A^_j.O“^wwlo ++ ■■•) +C’(as) +0{k‘^). (4.25) 


We already argued in refs. [65, 75] that an expansion of this kind up to and including 
0(5gK^) also captures the dominant N^LO QCD corrections, since the missing terms 
AloAjsisj^q, AloA^^j^q, and AJ^Aj^slo (represented by the dots) are small. It is there¬ 
fore justified to refer to the accuracy of our calculations as improved NNLO. Hence we 
expand all products in eq. (4.6) (and in all other observables) in 5s and k up to the afore¬ 
mentioned order, and neglect all higher terms. The observables are also expanded in the 
power-correction parameters Xi, 2 , Pi, fu^, fs up to linear terms. Higher powers as well as 
products of these parameters are dropped. 
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as{M^) = 0.1184 ±0.0007 

me = 0.51099892 MeV 

aeiM^) = 1/127.918 

m^ = 105.658369 MeV 

= sin^ Ow = 0.2312 

mr = 1.77699 GeV 

\VtlVtb/Vcb\^ = 0.9621 ± 0.0027 [94] 

mc{mc) = (1.275 ± 0.025) GeV 

= 130.5 ±11.6 [94] 

ml^ = (4.691 ± 0.037) GeV [95, 96] 

BRiB AceF)exp = 0.1051 ± 0.0013 [95] 

mt,poie = (173.5 ± 1.0) GeV 

Mz = 91.1876 GeV 

ms = 5.2794 GeV 

Mw = 80.385 GeV 

C = 0.574 ±0.019 [80] 

= 5t|5 GeV 

/io = 120+If GeV 

Xf = (0.12 ±0.02) GeV^ 

pi = (0.06 ± 0.06) GeV^ [97] 

Af = (-0.362 ± 0.067) GeV^ [95, 96] 

/o ± /. = (0 ± 0.2) GeV3 [60] 

/o - /, = (0 ± 0.04) GeV^ [60] 

ft = (0 ± 0.4) GeV3 [60] 


Table 1. Numerical inputs used in the phenomenological analysis. Unless specified otherwise, they 
are taken from PDG [98]. 


5 Phenomenological resnlts 

In this section we give the numerical results of our phenomenological analysis. We use the 
input parameters as given in table 1. For each variable we give the integral over bin 1 
(1 GeV^ < < 3.5 GeV^), bin 2 (3.5 GeV^ < < 6 GeV^), and the entire low-g^ region 

(1 GeV^ < q^ < 6 GeV^). In the high-g^ region we integrate over all q^ > 14.4 GeV^. The 
respective g^-interval is indicated by the argument of the observables. We give the numbers 
for electron and muon final state separately, and remind the reader that, depending on the 
channel and the experimental setup, our numbers have to be modified according to our 
Monte Garlo study in section 7. 

The quoted uncertainties are the parametric and perturbative ones only. Additional 
uncertainties from subleading non-perturbative corrections are not included. In particular, 
the 0{as{fib)-^(iCD/iTT‘c,b) non-perturbative corrections are estimated to be around ~ 5% in 
the low-g^ region. The individual error bars are obtained by varying the parameters in the 
range indicated in table 1, where we assume the errors on C and nic to be fully correlated. 
The total error is obtained by adding the individual ones in quadrature. By default we give 
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0[1,6] 

G [1,6] GeV2 

0(1,3.5] 

G [1,3.5] GeV^ 

'^0(1,3.5] '^0(1,3.51 

q^ 

0(3.5,61 

G [3.5, 6] GeV^ 

^0(3.5 6] ^ 0 ( 3 . 5,6] 

S[l,61 

6(1,6] 

0(1,61 

B[1,6] 

6(1,6] 

0(1,3.51 

6(1,6] 

6(1,6] 

0(3.5,61 

B 

100 

5.1 

5.1 

54.6 

3.7 

6.8 

45.4 

1.4 

3.1 

Ht 

19.5 

14.1 

72.5 

9.5 

8.8 

92.1 

10.0 

5.4 

53.6 

Ht 

80.0 

-8.7 

-10.9 

44.7 

-4.7 

-10.6 

35.3 

-4.0 

-11.3 

Ha 

-3.3 

1.4 

-43.6 

-7.2 

0.8 

-10.7 

4.0 

0.6 

16.2 


Table 2. Relative size of QED effects on 6 —>■ se+e“ at low-g^ (the muon case can be easily obtained 
by rescaling). All entries are given in percent. For each of the three bins the first two columns are 
the integrated observable and its QED correction normalized to the total low-g^ branching ratio, 
respectively Oj f^B and AO/ B). The third column is the relative size of the QED 

correction AO/ O). The sum of Ht and Hl does not exactly reproduce the branching 

ratio because in the latter we include finite bremsstrahlung and non-log enhanced QED corrections 
that are not available for the first two. 


two decimal digits. In case this leads to 0.00 we give the number up to the first significant 
digit. 

Before presenting our actual results, we would like to comment on the size of QED 
corrections. In table 2 the first two columns in each of the three sections are, respectively, 
the observable at NNLL and its QED correction expressed as a percentage of the branching 
ratio integrated in the whole low-g^ region (i3[l, 6]ee)- The third column is the relative size 
of the QED correction with respect to each NNLL observable. 

One can see immediately that the relative size of QED corrections to Ht is large, 
see third column in each section in table 2. Therefore, a few remarks on this observable 
are in order. It turns out that Ht is suppressed in the low-g^ region. To see this, let us 
look at the tree-level dependence of Ht and on the Wilson coefficients presented in 
eqs. (2.1) and (2.2). The phase space corresponding to \Cg\‘^ is suppressed in Ht compared 
to Ht, whereas that associated to \Ci\^ is enhanced. Surprisingly, this leads to a two-fold 
suppression of Ht- Eirst, there is an additional factor of 2s in the overall phase space w.r.t. 
Ht- Second, the factor \Cg + 2/sC^\^ is small in the low-g^ region, and even vanishes at 
the position of the zero of Ha- 

The QED corrections to the Hj, however, do not follow this pattern of suppression. 
In fact, from the inspection of the second columns in each section in table 2 we see that 
the absolute values of these corrections are natural in size and that all entries in these 
columns have the same order of magnitude. In the case of Ht the smallness of the NNLL 
QCD result implies that their relative size is anomalously large (see the third columns 
in table 2). However, we emphasize here that this does not indicate a breakdown of the 
perturbative series because the large relative size of QED corrections is almost entirely 
due to the suppression of the tree-level plus QCD contribution, and not due to a large 
absolute value of the QED corrections. To support our analytical findings, we investigated 
the situation in a Monte Carlo study (for details, see section 7) and find exactly the same 
pattern once we use EVTGEN and PHOTOS, see figures 12 and 13 in section 7.3. 
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We can even turn the argument around and regard the relative size of QED corrections 
in I-Lt as a virtue rather than a drawback, because it offers a good opportunity to be 
sensitive to QED corrections even without the pure QED observables and Hi defined 
in section 2. 

Finally, let us point out that similar large effects on Ha (or, equivalently, the forward- 
backward asymmetry) integrated in the whole low-q^ region are simply a due to the large 
cancellation between the integrated asymmetry in the two bins. This cancellation originates 
from the presence of a zero in the differential T-La spectrum and is not reproduced in the 
pattern of QED corrections. As we see in table 2, the latter imply a positive shift on %a 
in both bins. 

In the upper panel of figure 3 we show the differential distributions that we obtain 
for the various T-Li in the electron channel; dashed lines are obtained by switching QED 
corrections off. In the lower panel of figure 3 we show the log-enhanced QED correction 
itself, i.e. the difference between solid and dashed lines in the upper panel. 

5.1 Ht and Hl 

For the quantities Ht aad Tii we hnd theoretical uncertainties of 6 to 9%. In this sense 
the QED corrections listed in table 2 are really significant. 

3.5]ee =(2.91 ± O.lbscale ± 0.03mt ± 0.08c,mc =•= 0.02^4, 

± 0.003„, ± O.OlcKM ± 0.04 br,i) • 10"^ = (2.91 ± 0.19) • 10"^ , 

'Ht[ 3.5, 6]ee =(2.43 ± O.lGscale ± ± 0.08c,± 0.05^, 

± 0.01«, ± O.OlcKM ± 0.03 br,,) • 10-^ = (2.43 ± 0.20) • 10“^ , 

6]ee =(5.34 ± 0.33scale ± 0.07mt ± 0.16c,mc =•= O.OGmj, 

± O.Ola,, ± O.OlcKM ± 0.06 br,i) • lO”’^ = (5.34 ± 0.38) • 10“^ . (5.1) 

3.5]^^ =(2.09 ± O.lOscale ± 0.02mt ± 0.06c,me =t O.Olmt 

± 0.01„, ± O.OlcKM ± 0.03 br,,) • 10-^ = (2.09 ± 0.12) • 10“^ , 

77r[3.5, 6]^^ =(1.94 ± 0.13scale ± 0.03m4 ± 0.07c,me ± 0.05m, 

± 0.01„, ± O.OlcKM ± 0.02 br,,) • 10-^ = (1.94 ± 0.16) • 10"^ , 

77t[ 1, 6]^^ =(4.03 ± 0.23scale ± 0.06mt ± 0.12c,mc =t 0.06m, 

± 0.002„^ ± O.OlcKM ± 0.05 br,i) • 10“'^ = (4.03 ± 0.28) • lO""^ . (5.2) 


Tfc)!, 3.5]ee =(6.35 ± 0.23scale ± 0.08mt ± 0.22c,mc =t 0.08m, 

± 0.03a,, ± 0.02ckm ± 0.08br,i) • 10"^ = (6.35 ± 0.35) • 10"^ , 

'Hl[3.5, 6]ee =(4.97 ± 0.22scale ± 0.06m, ± 0.17c,me ± 0.04m, 

± 0.02a^ ± O.OlcKM ± 0.06br,J • 10-^ = (4.97 ± 0.29) • 10“^ , 
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9^(GeV^) 



9^(GeV^) 


Figure 3. Differential distributions for the various observables (upper panel) and their respective 
log-enhanced QED correction (lower panel) in units of 10“’^. Dashed lines correspond to switching 
off all QED corrections. The integrals under the curves reproduce the results presented in section 5 
for the electron channel. 


^l[ 1, 6]ee =(1.13 ± 0.04scale ± O.Olmt ± 0.04(7^^^ =•= O.Olm,;, 

± 0.01„^ ± O.OOScKM ± 0.01br,i) • 10“® = (1.13 ± 0.06) • 10“® . (5.3) 

3.5]^^ =(6.79 ± 0.23scaie ± 0.08^^ ± 0.23c,mc =•= 0.09^5 

± 0.03„, ± 0.02ckm ± 0.08br„) • 10“^ = (6.79 ± 0.36) • 10"^ , 

Hl[ 3.5, 6]^^ =(5.34 ± 0.23scaie ± 0.06^, ± 0.19c,± 0.04^, 

± 0.03„, ± O.OlcKM ± 0.07br,,) • 10-^ = (5.34 ± 0.32) • 10"^ , 

=(1-21 ± 0.04scale 4= O.Olmj zh 0.04c,mc =t O.Olm^ 

± 0.01„, ± 0.003CKM ± 0.02br,,) • 10-® = (1.21 ± 0.07) • 10“® . (5.4) 
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5.2 Ua 


For the zero-crossing Qq of T-La^ which is equivalent to the zero of the forward-backward 
asymmetry due to equation (4.3), we find 

{ql)ee =(3.46 ± O.lOscale ± 0.001^, ± 0.02c,± 0.06^, ± 0.02„J GeV^ 

=(3.46 ± 0.11) GeV^ , (5.5) 

=(3.58 ± O.lOscale ± 0.001^, ± 0.02c,± 0.06^, ± 0.02„J GeV^ 

=(3.58 ± 0.12) GeV^ . (5.6) 

We observe that the inclusive zero is in the same region as the semi-inclusive one obtained in 
the presence of a cut on mxs [70], but considerably lower than in the exclusive B -A- 
case [99]. The integrated T-La reads 

'HaI'^, 3.5]ee =(-1.03 ± 0.04scale ± O.Ol^t ± 0.02c,mc =•= 0.02^^ 

± 0.01a, ± 0.003ckm ± 0.01br,i) • 10 "^ = (-1.03 ± 0.05) • 10 "^ , 

-^^[3.5,6]ee =(+0.73 ± O.llscale ± 0.01^, ± 0.04c,m, ± 0.05^, 

± 0.02a, ± 0.002 ckm ± O.OIbr,,) • 10"^ = (+0.73 ± 0.12) • 10"^ , 

77a[ 1, 6]ee =(—0.29 ± 0.14scale ± 0.002^4 ± 0.02c,mc ± 0.06mj 

± 0.03a, ± O.OOlcKM ± 0.004 br,i) • lO”"^ = (-0.29 ± 0.16) • 10“'^ . 

(5.7) 

77a[1, 3.5]^^ =(-1.10 ± 0.03scale ± O.Olmt ± 0.02c,m, ± 0.02^6 

± 0.01a, ± 0.003ckm ± O.OIbr,,) • 10 "^ = (-1.10 ± 0.05) • 10 "^ , 

77a[3.5, 6 ];,;, =(+0.67 ± O.llscale ± O.OU, ± 0.04c,m, ± 0.05^, 

± 0.02a, ± 0.002 ckm ± O.OIbr,,) • 10“^ = (+0.67 ± 0.12) • 10“^ , 

77a[1, 6 ]^^ =(—0.42 ± 0.14scaie =•= 0.003^* ± 0.01c,mc =•= 0.06^^ 

± 0.03a, ± O.OOlcKM ± 0.01 br,i) • lO”’^ = (-0.42 ± 0.16) • 10“^ . 

(5.8) 

As far as the total error is concerned, the single bins are much better behaved than the 
entire low-g^ region. This is due to the large cancellation of the central values of bin 1 and 
bin 2, which is owed to the presence of the zero. The value of the latter happens to be 
almost exactly at the position where we subdivide the low-q^ region into bin 1 and bin 2. 

5.3 773 and 774 

For the observables 773 aad 774, sensitive to QED corrections, we find 

773[1, 3.5]ee =(4.04 ± 0.64scaie ± 0.04^, ± 0.13c,™, ± 0.10™, 

± 0.03a, ± O.OlcKM ± 0.05 br,,) • 10-9 = (4.04 ± 0.67) • lO'^ , 
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?^3[3.5, 6]ee =(4.88 ± O.SOscale ± 0.05^, ± 0.16c,m, ± 0.07^, 

± 0.02«, ± O.OlcKM ± 0.06br,i) • 10-9 = (4.88 ± 0.54) • lO"® , 

'^3[1, 6]ee =(8.92 ± 1.14scale ± O.lOm^ ± 0.30c,me 4= O.lGmj 

± 0.06«, ± 0.03CKM ± O.IIbr,,) • 10-9 = (8.92 ± 1.20) • 10-9 . (5.9) 

3.5]^/^ =(1.68 ± 0.26scale 4= 0.02mt 4= 0.06c,mc 4= 0.04mfc 

± 0.01a, ± 0.005CKM ± 0.02br„) • 10-9 = (1.68 ± 0.27) • 10-9 , 

773[3.5,6]^/i =(2.03 ± 0.21scaie 4= 0.02mt 4= 0.07c,mc 4= 0.03mb 

± 0.01a, ± 0.006ckm ± 0.03br,i) • 10-9 = (2.03 ± 0.22) • 10-9 , 

'^3[1, =(3.71 ± 0.47scale 4= 0.04^* 4z 0.12c,me 4= 0.06mfc 

± 0.02a, ± O.OlcKM ± 0.05br„) • 10-9 = (3.71 ± 0.50) • 10-9 . 

(5.10) 

774[1, 3.5]ee =(6.23 ± 0.55scale ± 0.07m, ± 0.21c,me ± 0.01m, 

± 0.02a, ± 0.02CKM ± 0.08br,J • 10-9 = (6.23 ± 0.60) • 10-9 , 

'H4[3.5, 6]ee =(2.19 ± 0.16scale 4= 0.03mt 4= 0.07c,mc =•= 0.02m, 

± 0.006a, ± O.OOGckm ± 0.03br,i) • 10-9 = (2.19 ± 0.18) • 10-9 , 

'H4[1, 6]ee =(8.41 ± 0.71scale 4= O.lOm, 4= 0.28c,me =•= 0.02m, 

± 0.02a, ± 0.02CKM ± O.IObr,,) • 10-9 = (8.41 ± 0.78) • 10-9 . 

(5.11) 

'H4[1, 3.5]/,^ =(2.59 ± 0.23scale 4= 0.03mt 4z 0.09c,me 4= 0.006m, 

± 0.007a, ± 0.007 ckm ± 0.03br„) • 10-9 = (2.59 ± 0.25) • 10-9 , 

'H4[3.5, =(0.91 ± 0.07scale 4= O.Olmt 4z 0.03c,me 4= 0.008m, 

± 0.002a, ± 0.003 ckm ± O.OIbr,,) • 10-9 = (0.91 ± 0.075) • 10-9 , 

774 ( 1 , =(3.50 ± 0.29scale 4= 0.04mt 4z 0.12c,me 4= 0.01m, 

± 0.01a, ± O.OlcKM ± 0.04br„) • 10-9 = (3.50 ± 0.32) • 10-9 . 

(5.12) 

5.4 Branching ratio, low-g^ region 

The decay width is simply given by the sum of Ht and Hi and hence can in principle be 
derived by the numbers given in the previous subsections. However, we give the numbers 
explicitly here, for two reasons. First, the branching ratio is an important quantity, also 
experimentally. Second, there are two more contributions which are available only for the 
branching ratio, but not for "Hr or 77 l individually. These are the finite bremsstrahlung 
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contributions from [47] and the non-log enhanced terms of Wgg™^(s). Both give only a small 
correction, but we include them for the sake of completeness. This yields 


B[l, 3.5]ee =(9.26 ± 0.34scale ± O.IU, ± 0.30c,m, ± 0.10^, 

± 0.02„^ ± 0.03ckm ± O.IIbr,,) • 10“^ = (9.26 ± 0.49) • lO'^ , 

B[3.5, 6]ee =(7.44 ± 0.37scale ± 0.10^, ± 0.26c,± 0.08^, 

± 0.03„, ± 0.02ckm ± 0.09br,,) • 10“^ = (7.44 ± 0.48) • lO'^ , 

B[l, 6]ee =(1.67 ± 0.07scale ± 0.02^, ± 0.06c,m, ± 0.02^, 

± 0.01a,, ± O.OOScKM ± 0.02br,i) • 10"® = (1.67 ± 0.10) • 10"® . (5.13) 

B[l, 3.5]^^ =(8.88 ± 0.31scaie ± 0.11^, ± 0.29c,m, ± 0.10^, 

± 0.02a^ ± 0.02ckm ± O.IIbr,,) • 10"^ = (8.88 ± 0.46) • 10"^ , 

B[3.5,6]^^ =(7.31 ± 0.36scaie ± 0.09^, ± 0.25c,m, ± 0.09^, 

± 0.03a, ± 0.02ckm ± 0.09br„) • 10"^ = (7.31 ± 0.47) • 10"^ , 

B[l, 6]^^ =(1.62 ± 0.07scale ± 0.02mt ± 0.05c,mc =t 0.02mj, 

± 0.01a, ± O.OOScKM ± 0.02br,i) • 10"® = (1.62 ± 0.09) • 10"® . (5.14) 

The values are about 2% larger compared to our previous analysis [75]. This is due to 
updated input parameters and the inclusion of the Kriiger-Sehgal corrections [71, 72]. 

5.5 Branching ratio, high-g^ region 

The branching ratio in the high-g^ region suffers from large uncertainties stemming from 

2 3 

hadronic input parameters in the l/m^’ power-corrections, which results in total error 
bars of 0(30%), 


B[> 14.4]ee =(2.20 ± 0.30scale ± 0.03mt ± 0.06c,mc =•= 0.16^^ ± 0.003a, ± O.OlcKM ± 0.03br,i 
± 0.12a 2 ± 0.48pi ± 0.36/, ± 0.05/J • 10“^ 

=(2.20 ±0.70) • 10“^ , 

B[> 14.4]pp =(2.53 ± 0.29scaie ± 0.03,„, ± 0.07c,m, ± 0.18,„, ± 0.003a, ± O.OIckm ± 0.03br„ 

± 0.12a 2 ± 0.48pi ± 0.36/, ± 0.05/J • 10“^ 

=(2.53 ± 0.70) • 10" J (5.15) 

Comparing these results to earlier analyses on the high-( 7 ^ branching ratio shows that 
our numbers are considerably lower than the ones in [52, 53]. In the following, we show 
that this is the result of several effects which all give corrections in the same direction. 
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Once we turn to the prescriptions given in [52, 53] we reproduce their results, as can be 
seen below. 

We first perform the comparison to Greub et ah [53]. We start with the above num¬ 
bers and first switch off the In-enhanced QED corrections, which also removes the 
difference between the muon and the electron channel, and yields 2.74 (all numbers that 
follow are in units of 10“^). Next, we turn off the finite bremsstrahlung contributions, 
which is only a minor effect and does not change the digits given before. Taking out the 
Kriiger-Sehgal corrections, on the other hand, is a rather large effect in the high-g^ region 
and results in 3.05. We also have to remove the 1/m^’^ and l/m^ non-factorisable power- 
corrections which further increases the result to 3.36. Switching furthermore off those QED 
corrections which are not ln(m^/m£)-enhanced, we get 3.56. This shift is rather large, but 
we remind the reader that some of these terms are enhanced by sin^ Gw). Chang¬ 

ing from four- to two-loop running for has again only a minor impact and gives 3.55. 
We now switch off the change in renormalisation scheme for the quark masses, i.e. we use 
the pole mass for charm and bottom. Eurthermore, we use the input parameters from [53]. 
Both effects taken together give 3.68. We now take into account that the integration in¬ 
terval in [53] is given in the variable s = /m^. Hence a change in the value for mb 

results in the modified lower integration limit = 13.824 GeV^. This effect must not 
be underestimated because it brings the branching ratio up to 4.36. We now turn to the 
normalisation prescription given in [53], which instead of the factor C from eq. (4.5) and 
the perturbative expansion of r(6 —)• uev) makes direct use of the perturbative expan¬ 
sion of r(6 — ceu), including charm-mass dependent phase-space factors and radiative 
corrections. This increases the branching ratio further to 4.57. Einally, we divide by the 
experimentally measured semileptonic 5 —)• c branching ratio (see table 1) and get 43, which 
is precisely the value of i?high, pert in eq. (48) of [53] . 

The comparison to Ghinculov et al. in [52] proceeds along the same lines. The differ- 
ences to the analysis by Greub et al. are the Kriiger-Sehgal corrections and the , 

1/m^ power corrections, both are taken into account in [52]. Moreover, different input pa¬ 
rameters are used and the lower integration limit is formulated in rather than in s. To 
quantify these effects, we first switch off again ln(m^/m£)-enhanced QED corrections and 
finite bremsstrahlung effects first and end up with 2.74. We then also remove those QED 
corrections that are not enhanced by ln(m^/m^), which gives 2.93. Ghanging from four- 
to two-loop running for Ug is again only a small effect and gives 2.92. The biggest effect 
comes from the change of input parameters and the removal of the renormalisation-scheme 
conversion for the quark masses, i.e. we now use the pole mass for charm and bottom. 
These two effects taken together result in 3.89. Einally, we switch to the normalisation 
that is used in [52] and get 4.02. This number coincides within a fraction of a percent with 
the value 4.04 from eq. (6.36) in [52]. The obtained level of accuracy shall be sufficient for 
the present check. 

5.6 The ratio TZ{sq) 

77(14.4)ee =(2.25 ± 0.12scale ± 0.03^^ ± 0.02c,nic =t O.Olm^ ± O.Ola^ ± 0.20cKM 
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± 0.02a, ± O.IV ± 0.08^0+^^ ± 0.12^o_^J • lO'^ 

= (2.25 ±0.31) • 10"^ , 

'^(14.4)^^ =(2.62 ± O.ODscale ± 0.03mt ± 0.01c, me =•= O.Olm^ ± O.OIq,^ ± 0.23ckm 
± 0.0002a, ± 0.09p, ± 0.04^0+^^ ± 0.12^o_^J • lO’^ 

= (2.62 ± 0.30) • 10“^ . (5.16) 

We clearly see a reduction of the total error bars from 0(30%) in the high-q^ branching 
ratio to 14% and 11% in the electron and muon channel of TZ{sq), respectively. Besides the 
uncertainties due to power corrections, also the scale uncertainty gets significantly reduced. 
The largest source of error are CKM elements (notably Vub)- 


6 New Physics sensitivities 


In this section we present the constraints on the most relevant Wilson coefficients (Cg and 
Cio) that we obtain using the current experimental results, and investigate the reach of 
Belle II with an expected final integrated luminosity of 50 ab~^. 

Previous model-independent new physics analyses [76, 100-102], as well as studies in 
specific models such as minimal-flavour-violation [103-105], two-Higgs doublet models [106, 
107], and supersymmetry [101, 108-114] can be found in the literature. 

The weighted averages for the low- and high-g^ branching fractions have been presented 
in eq. (1.2). Here we need the results on the individual channels: 


BiB ^ )“P 


’ (1.493 ± 0.504l°J^|) X 10-6 (Belle, U) 
(l.93l[];^^;[];?^ ± 0.18) x lO"® (BaBar, ee) 
(0.66l[];^^;[]|^ ± 0.18) X lO”® (BaBar, ^u/r) 
± 0.18) X 10-6 (BaBar, U ), 


( 6 . 1 ) 


B{B 



(0.418 ±0.117+|[;[J^^) X 10-6 
(0.56l6;J™3 ± 0.00) X 10-6 
(0.6016;6116;06 ± o.OO) x 10-6 
_ (0.57l6;J610;03 ± 0.00) X 10-6 


(Belle, U) 
(BaBar, ee) 
(BaBar, fifj, ) 
(BaBar, U ). 


( 6 . 2 ) 


In each result, the first error is statistical, the second systematics and the third model- 
depedent systematics which is included in case of Belle in the second error. Note that the 
high-g^ region chosen by BaBar and Belle have a slightly different minimum (14.4 and 
14.2 GeV^ for Belle and BaBar, respectively). 

In ref. [38] Belle presented a measurement of the normalized forward-backward asym¬ 
metry defined in eq. (2.10) in the low- and high-g^ regions. The binning chosen to present 
the measurement (binl = [0.2,4.3] GeV^ and bin2 = [4.3,7.3(8.1)] GeV^ for electrons 
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(muons)) differs from the one proposed in this work. In particular, the larger integration 
end-point in the second bin includes a region of the spectrum where sizable interference 
from the tail of the J/i!) is present. From ref. [38] we read: 


-exp,- , . 0.34 ±0.24 ±0.02 Belle, binl 

4b B ^ x,tr) = { 6 

[0.04 ±0.31 ±0.05 Belle, bin2. 

In order to preserve the cancellation of systematic uncertainties. Belle averaged the nor¬ 
malized asymmetries in the electron and muon channels; i.e. 

Ayb{B ^ Xstr) = ^ Xse+e-)+ABBiB ^ X,4^-)) /2 . (6.4) 


We integrated our differential spectra in the above bins in order to investigate the impact 
that this measurement has on the Wilson coefficients, but we caution the reader that the 
uncertainties we quote could be underestimated. We find: 


Abb{B ^ Xs£+r) 


-0.0773 ± 0.0057 binl 
±0.049 ± 0.018 bin2. 


(6.5) 


We define the following ratios of high-scale Wilson coefficients (see [75] for the precise 
definitions of the Wilson coefficients). 




and 


7 ? 9,10 


^9,10 (t^o) 


( 6 . 6 ) 


The numerical formulas for all observables in terms of the ratios Ri can be found in 
appendix B. We assume that the relative theoretical uncertainty on a given observable 
{60/0) is mostly independent of the precise values of Wilson coefficients and that it can 
be extracted from the SM predictions presented in section 5. 

We present the bounds on the ratios Rg and Rio under the assumption of no new 
physics contributions to the magnetic and chromo-magnetic dipole operators {Rt^s = 1) 
in figure 4 (similar analyses were done, e.g., in [76, 100]). The contours are the 95% C.L. 
regions allowed by the experimental results in eqs. (6.1), (6.2) and (1.2); two sigma theo¬ 
retical uncertainties are added linearly. In each plot we show the impact of the branching 
ratio measurement in the low-g^ (red regions) and high-f;^ (green regions) and their overlap 
(black regions). The SM corresponds to the point [i? 9 ,i?io] = [1,1]. As we discuss below, 
the small yellow contours correspond to the Belle II estimated reach, assuming that the 
observed central values coincide with our predictions. The top left, top right and lower 
plot consider the B —)> Xse~^e~, B —)• XsfJ,'^^~ and B —)> Xsi'^i~ cases, respectively. In 
the lower plot in figure 4 we include also the 95% C.L. bounds from the Belle measurement 
of the normalized forward-backward asymmetry given in eq. (6.3); the region outside the 
dashed and dotted parabola shaped regions are allowed by the measurement in binl and 
bin2, respectively. The resulting picture is in overall agreement with the SM expectations at 
the 95% C.L.; though we should note that at the one sigma level there are some statistically 
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Figure 4. Constraints on the high-scale Wilson coefficient ratios [i? 9 ,i?io] that we obtain at 95% 
C.L. from the present BaBar and Belle experimental branching ratios measurements. In the upper 
left (upper right, lower) plot we show the constraints obtained from the measured branching ratios in 
the low-g^ and high-g^ region in the electron (muon, electron plus muon) channel. The red and green 
regions correspond to the low- and high-g^ regions, respectively. The black region is the overlap of 
these two constraints. The dot is the SM expectation ([i? 9 , i?io] = [1,1]). The yellow contour is the 
Belle II reach (see figures 5-7). The region outside the dashed (dotted) parabola shaped regions are 
allowed by the Belle measurement of the normalized forward-backward asymmetry in binl (bin2). 


insignificant tensions driven by a disagreement between low- and high-g^ measurements in 
the muon channel. 

In order to study the expected Belle II reach, we estimate the statistical uncertainties 
on the various observables using the squared weight method detailed in ref. [115]. Let us 
consider the following differential quantity: 

_ c up^od 

dsdz Ftot dsdz 


where the £ is the integrated luminosity, CTprod is the production cross section for e+e —>• 
BB at the B-factories’ center of mass energy, Btot is the total B decay width and d'^T/dsdz 
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[1,3.5] 

[3.5,6] 

[1,6] 

> 14.4 

B 

3.7 % 

4.0 % 

3.0 % 

4.1% 

Bt 

24 % 

21 % 

16 % 

- 

Bl 

5.8 % 

6.8 % 

4.6 % 

- 

Ba 

37 % 

44 % 

200 % 

- 

B 3 

240 % 

180 % 

150 % 

- 

Bi 

140 % 

360 % 

140 % 

- 


Table 3. Statistical uncertainties that we expect at Belle II with 50 ab ^ of integrated luminosity. 
The first row gives the considered bin in GeV^. 


is the double differential B —)> decay rate. The number of events that we expect 

to observe in a certain range of s and z is 


3Xp - 

f 'll • 

J dsdz 

(6.8) 

<^Wexp = 

's/ 

(6.9) 


where (5A4xp is the expected statistical error. If instead of considering simple slices of the 
integration region we utilize a weight function IT[.s, z] to define an observable (that cannot 
be anymore interpreted in terms of “number of events”), the above equations generalize to 


Oexp = 
•JOexp = 



W\s, z\ ds dz , 



1 

IT[s, z]^ ds dz 


( 6 . 10 ) 

( 6 . 11 ) 


Note that eq. (6.11) reproduces the correct uncertainties for the simple case in which the 
weight is a product of theta functions (i.e. the integral is restricted to a certain region of 
phase space) and that the relative uncertainty (5C>exp/Oexp is invariant under rescaling of 
the weight function. 

In order to asses expected uncertainties on observables corresponding to the weights 
given in eq. (2.7) we start from the double differential rate given in eq. (1.4) and use the 
expressions for the Hj in eqs. (2.1)-(2.3) and use some reference value for the Wilson 
coefficients. Next we fix the normalization £ fprod/htot in such a way to reproduce the 
~ 25% statistical uncertainty that BaBar obtains with an integrated luminosity £current = 
0.4242 ab“^ [36]. Finally we rescale the normalization by the factor £future/23current where 
23future = 50 ab“^ is the Belle II expected final integrated luminosity. 

This procedure produces acceptable error estimates for Bt, Bl and Ba, while fails 
for Bs and B 4 . The reason is that the integral in eq. (6.10) vanishes when integrated the 
simple NLO formula given in eq. (1.4) against the weights W^a- We bypass this problem 
by extracting Oexp from the exact results presented in section 5 and using eq. (6.11) to 
calculate the error (in fact the weights (IF 3 ^ 4 )^ do not annihilate the NLO differential 
width). 
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Figure 5. Constraints on [Rg, i?io] that we expect at 95% C.L. from Belle II measurements of the 
branching ratio in the low-g^ (left plot) and high-g^ (right plot) regions with 50 ab“^ of integrated 
luminosity. For the low-g^ case, the solid and dashed contours correspond to the branching ratio 
restricted to the low ([1,3.5] GeV^) and high ([3.5,6] GeV^) bin, respectively. The hashed region is 
the overlap of the expected constraints from these two bins. The shaded region is the constraint we 
obtain by considering the branching ratio integrated in the whole low-g^ region. The black dot is 
the SM expectation. The solid red area is the overlap of all constraints we consider (it corresponds 
to the yellow contour in figure 4). 




Rg Rg 

Figure 6. Constraints on [i? 9 ,i?io] that we expect at 95% C.L. from Belle II measurements of 
T-Lt (left plot) and Rl (right plot) in the low-g^ region with 50 ab“^ of integrated luminosity. See 
figure 5 for further details. 


Following the discussion summarized in section 5 of ref. [17], we add a flat 2% systemat- 
ics to the projected statistical errors obtained with the squared weight method and obtain 
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Figure 7. Constraints on [Rg, i?io] that we expect at 95% C.L. from Belle II measurements of Ha 
(upper plot), H 3 (lower left plot) and H 4 (lower right plot) in the \ow-q^ region with 50 ab“^ of 
integrated luminosity. See figure 5 for further details. 


the low-g^ uncertainties collected in table 3. The expected uncertainty on the high-g^ 
branching ratio is taken directly from ref. [17]; in fact, near the end-point of the spectrum 
our method fails to take into account the improvement in the signal-to-background ratio. 

In figures 5, 6 and 7 we show the expected impact of Belle II measurements on the 
various observables we consider in the [i? 9 ,iiio] plane. Each contour is drawn at 95% 
C.L. by combining linearly theoretical and experimental uncertainties. In the scenario we 
consider the strongest bounds on the Wilson coefficients are driven by measurements of 
the low-q^ branching ratio and of Ha and Ht in the two bins. The latter statement is 
driven by the assumption that the future experimental central values will coincide with the 
respective SM expectations. If deviations are seen, all observables become crucial to pin 
down the structure of new physics. 
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Figure 8. Pictorial descriptions of the theoretical definition (left) of B ^ and of the 

experimental recoil technique. 


7 On the connection between theory and experiments 
7.1 Various experimental settings 

Here we discuss how to compare integrated low and high-g^ observables, calculated with 
the inclusion of QED corrections, to quantities measured by BaBar, Belle and also by 
the future Belle II experiment. As we explain below, we find that our results can be 
directly compared to integrated observables measured at BaBar, Belle, and Belle II with 
the exception of the di-electron case at BaBar. In the latter case we have to increase our 
predictions for the integrated branching ratio in the low (high) region by 1.65% (6.8%), 
see eqs. (7.1) and (7.2). 

From the theoretical standpoint the Xg system, in the inclusive Xat'^t~ final state, 
contains all the electromagnetic radiation produced in the hard interaction, see the diagram 
on the left in figure 8. From the experimental point of view there are two distinct techniques 
to measure the inclusive B —?■ rate: the recoil and sum-over-exclusive methods. 

In the recoil technique, whose luminosity requirement makes it viable only at super flavor 
factories, one of the B mesons produced in the e'’'e“ hard interaction is tagged using a 
semileptonic or hadronic decay and the final state is identified by the two leptons only, see 
the diagram on the right in figure 8. In the sum-over-exclusive method, the recoling heavy 
meson is not looked at and the decaying B is fully reconstructed in final states with a 
and up to four pions. The fully inclusive rate is then reconstructed using JFTSFT [116]. 

The comparison between the measured branching ratio (BR) and the results of our 
inclusive calculations depends critically on the definition of (f‘. If no photons are included 
in the definition of the di-lepton invariant mass (i.e. (f' = {pi+ + Pe-)'^) our results can be 
used directly in the comparison with experiments. This is the case for the di-muon channel 
at both experiments [117, 118] and for the di-electron channel at Belle [117]. This will be 
exactly the case in a fully inclusive analysis using the recoil technique at Belle II. However, 
at BaBar photons that belong to a B —)• Xse^e~ event and that are emitted in a cone of 
35 mrad angular opening around either final state electron are included in the calculation 
of the [118]. 
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In order to calculate the shift that the latter definition has on to the inclusive theory 
prediction we generate inclusive B —)• events using EVTGEN [119], hadronize them 

with JETSET and include electromagnetic radiation with PHOTOS [120, 121], Following 
the BaBar and Belle procedure we build a fully inclusive sample in the whole and 
mXs phase space by fully inclusive events (parton level supplemented by a Fermi Motion 
Model [56]) for mx^ > 1-1 GeV with exclusive B —?■ events (to describe the 

low mxs region). Using this large event sample we were able to calculate the impact of 
including photons emitted in a 35 mrad cone around either electron in the q'^ calculation. 
We find: 


~lsi 


g=Pe++Pe-+P7coll 

owl 

ee lq=p +-\-p 


- 1 = 1.65% 






- 1 = 6.f 


-l9=Pe++Pe- 


(7.1) 


(7.2) 


where the suffixes q = Pg+ + Pg- a-iid Q = Pe+ + Pe- + 7*70011 refer to quantities we calculate 
and observables measured at BaBar, respectively. 


7.2 Validation 

The results presented in the previous subsection depend crucially on the reliability of using 
PHOTOS to model photon radiation in 6 —?■ decays. In this subsection we perform 

several checks to validate this approach; in particular we show that PHOTOS can be used 
to reproduce (to a good enough extent) the effects of QED radiation that we calculate 
analytically. 

As discussed above, we generate inclusive B —)• Xs£'^i~ events using EVTGEN, 
hadronize them with JETSET and include electromagnetic radiation with PHOTOS. In 
order to obtain a fully inclusive event set we combine K, K* and Xs{mxs > 1-1 GeV) 
samples. The mxs and q^ spectra that we obtain are presented in figure 9. The relative 
weights of the K and K* samples with respect to the inclusive {mx^ >1.1 GeV) one have 
to be provided externally. The actual weights we adopt are extracted from experimental 
results for the exclusive and inclusive modes and their precise values do not impact much 
the shape of the q^ spectrum. In fact, as we can see in the plot on the right of figure 9 
only the very high di-lepton invariant mass region, q"^ > 17 GeV^, is affected. 

A point that is important to mention is that PHOTOS generates events with large 
photon multiplicity while analytic calculations are confined to a single photon emission. 
Obviously the vast majority of photons emitted are soft and/or collinear to the final state 
leptons; moreover, only relatively high energy collinear photons can impact the shape of 
the spectrum. 

In the left panel of figure 10 we show the photon multiplicity we observe in the gener¬ 
ated events. The shaded area corresponds to events for which the most energetic photon 
has < 30 MeV and that, at the experimental level, are identified as purely hadronic 
B —)■ Xsi'^£~. As expected there is a very large multiplicity of soft photons. We find that 
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Figure 9. mx^ and spectra that we obtain in a -B —?► Xsi~^£ sample generated combining the 
exclusive B —>■ K^*^£~^£~ modes with a pure inclusive calculation for mx^ > 1-1 GeV. 


-S2 

s; 


60 

50 

40 

30 

20 

10 

® 0 1 2 3 4 5 

Number of Photons 


I I All photons 
I I Photons with E <30 MeV 




Figure 10. Left: Distribution of events with < 5. Right: distribution of the most energetic and 
second most energetic photons. 


only 17% of all events (this is the integral of the purple unshaded region) correspond to 
final states with at least one photon with energy larger than 30 MeV. These photons are 
resolved experimentally and need to be included in the hadronic (V^) or leptonic 
system least the event is rejected (cf. also the last paragraph of this subsection). 

In the right panel of figure 10 we show the distribution of the most and second most 
energetic photon. The integral of the upper (lower) curve over a photon energy range 
[E^i,Ey 2 ] yields the percentage of events in which the most (second-most) energetic photon 
has energy in that interval. The fraction of events with at least one (two) photons is 58% 
(23%), is given by the integral of these curves and can also be easily read off from the 
left panel of figure 10. Since the impact of including certain collinear photons in the 
definition of the is more pronounced for more energetic photons, we see that these 
effects are completely described by a single photon emission: the analytic calculation of 
QED radiation is, therefore, completely adequate to discuss this phenomenon. 

Finally, in order to verify whether PHOTOS correctly models photon radiation in this 
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Figure 11. Effect of the inclusion of electromagnetic radiation calculated using EVTGEN + 
PHOTOS (left) and using analytical methods (right). 


decay, we need to compare spectra calculated with and without the inclusion of QED 
radiation. Therefore, we generated a second set of events in which we switched PHOTOS 
off. The result of this analysis is presented in figure 11. In the left and right panels 
we show the Monte Carlo study and the result of our analytical calculation, respectively. 
Numerically, the relative shifts that we obtain for the branching ratio in the low and high-g^ 
regions are (in round brackets we present the analytical results): 


6BK{B ^ Xs^l^li~) 

6 BR{B X^e+e-) 


j +1.5%(+2.0%) low 
4.4%(—6.8%) high 

+3.6%(+5.2%) low 
-12.9%(-17.6%) highg2 


(7.3) 

(7.4) 


Given the differences in the techniques used, the agreement is remarkable. We conclude 
that the PHOTOS description of electromagnetic radiation is sufficiently close to the exact 
calculation to be used to reliably calculate the shifts we presented in eqs. (7.1) and (7.2). 

Before concluding this subsection, we would like to stress that validating the use of 
PHOTOS is important in its own right because experiments use it to estimate the impact 
of missing photons on their efficiencies. Legitimate B —>• Xs(l.~^l~ events might be rejected 
because of two possible reasons. First, if a large number of soft photons {E^ < 30 MeV 
and 20 MeV for BaBar and Belle, respectively) is present, they might push the event out 
of the mES^ and XE acceptance windows (see, for instance, refs. [35, 37] for a definition 
of these kinematical quantities). Second, if a photon with energy larger than 30 (20) MeV 
is not identified, most likely the event is discarded because the total momentum fails to 
reconstruct a decaying B meson. The latter effect can be quite substantial because, as we 
discussed above, about 17% (18%) of all B —)• Xs£^i~ events have at least one photon 

^Belle names this quantity mtc- 
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Figure 12. q^-dependence of Ht, 'Hl and branching ratio {Ht + Hl) that we extract from a 
B —>■ Xs£~^£~ sample generated combining the exclusive B —>■ K^*'>£'^£~ modes with a pure inclusive 
calculation for mx, >1.1 GeV. The dotted lines are obtained by switching off QED radiation. 


with energy larger than 30 MeV (20 MeV). The fraction of events that is lost to these two 
mechanisms is taken into account, in the calculation of the efficiencies, using PHOTOS. 


7.3 Monte Carlo estimate of QED corrections to %t and 


The results presented in section 5.1 indicate that the relative size of QED corrections to 
1 -Lt are about an order of magnitude larger than the corresponding corrections to 1 -Ll and 
to the branching ratio. In this section we show that this result is actually reproduced in 
our Monte Carlo study. As a first step we plot in figure 12 the spectra for Ht, 
and the branching ratio with (solid lines) and without (dotted lines) the inclusion of QED 
radiation. 

Note that the absolute size of QED effects on T-Lt, T~Ll and 'Ht + 'Hl is very similar and 
natural in size; in particular, a small positive net contribution to the integrated branching 
ratio in the low-gr^ region is the sum of a small negative shift on TIl and a slightly larger 
positive shift on T-Lt- We plot the actual QED corrections to the three observables in 
figure 13. 

From inspection of the left plot in figure 12 we see that, in the low-q^ region T-Lt is 
much smaller than Tii. We can understand the origin of this effect by looking at the ratio 
'Ht/'Hl at leading order: 


Hl ^Cfo + {Cg + 2C^y 


(7.5) 


The suppression comes from the small 2s < 1 factor and from the accidental strong can¬ 
cellation between Cg and 2 CjIs at low s (in fact, the combination Cg + 2 CjIs vanishes for 
s ~ 0.15). In the Standard Model Cy is negative; if its sign was reversed we would obtain 
Cg + 207/3 > Cg + 2 C 7 and the integrated T-Lt and Hl observables at low-( 7 ^ would assume 
very similar values. 

In table 4 we present the results we obtain by integrating the Monte Carlo generated 
b —)■ s££ histograms. For each bin ([si,S 2 ]) and for each observable O {Ht + 
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Figure 13. Differential distributions of the QED corrections to Ht, 'Hl and branching ratio 
{I-Lt +'Hl) that we obtain in a B ^ Xs£'^£~ sample generated using EVTGEN and PHOTOS and 
combining the exclusive B —>■ modes with a pure inclusive calculation for mx^ > 1-1 GeV. 
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1.0 

2.5 

Ht 
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5.0 
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3.0 

36.0 

Hl 

81.0 

-4.5 

-5.5 
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-2.5 

-5.0 

35.0 

-2.0 
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Table 4. Relative size of QED effects at low-g^ that we extract from our Monte Garlo b —?► 
se~^e~ sample (All entries are given in percent). For each of the three bins the two columns are 
the integrated observable and its QED correction normalized to the total low-q^ branching ratio 
O/ B and AO/ B). The third column is the relative size of the QED correction 


and T-Ll) we show the total integrated observable 0 / + 'Hl)), the total inte¬ 
grated QED effect AO/ +'Hl)) and the relative size of the QED correction 

AO/ O). We see that the absolute size of QED corrections is very similar amongst 
the three observables (with the effect on Ht being only slightly larger) and that the sup¬ 
pression of Ht with respect to Hl is responsible for very large relative effects in the 30-50% 
range. 

Einally we must point out that the numerical estimates presented in table 4 are affected 
by sizable uncertainties that are hard to quantify and that only the analytical results 
presented in table 2 should be utilized. The Monte Carlo study was nevertheless extremely 
valuable to build confidence in our study. 


- 39 - 






















8 Conclusion 


The inclusive decay B —)• is one of the most important modes in the indirect 

search for new physics via quark flavour observables. It is theoretically clean, while the 
exclusive mode is affected by unknown power corrections. Thus, besides allowing for a 
nontrivial check of the recent LHCb data on the exclusive mode, it contains complementary 
information both in Standard Model predictions and in pinning down new physics. It is 
therefore a precious channel to be measured at Belle II, and might be accessible even at 
LHCb. 

In the present article we perform a complete angular analysis of the inclusive decay 
B —)■ Xsi~^£~ by taking into account all perturbative and power corrections that are avail¬ 
able to date. We confirm the findings of ref. [76] that a separation of the double differential 
decay width into three observables as well as subdivision of the low-( 7 ^ region 

into two bins (see also [65]), provides significantly more information than the branching 
ratio or forward-backward asymmetry in the entire low-g^ region alone. 

We compute logarithmically enhanced QED corrections to these observables and find 
that they do not obey the simple second-order polynomial in z = cos(O) exhibited by the 
double differential decay width in the absence of QED corrections. We therefore propose 
to project out using weight functions, and argue that the Legendre polynomials 

Pn{z) are the optimal choice for the latter. Besides reproducing Hxiq^) and HL{q^) in 
the absence of QED radiation, they allow to construct observables 773 , 4 ( 0 '^) (eq. (2.6)) 
that vanish if only QCD corrections are taken into account, and are therefore particular 
sensitive to QED effects. In view of the benefits of the Legendre weight functions we 
urgently recommend the experiments to use the weights (2.6) to extract single-differential 
distributions, and to refrain from attempting polynomial fits to the data. 

The absolute values of the QED effects that we compute are natural in size. However, 
due to the phase-space and Wilson coefficient suppression of Hxiq'^) the relative size of the 
QED corrections is large in this observable. We argue carefully that this does clearly not 
indicate a breakdown of perturbation theory. On the contrary, we can benefit from the 
fact that QED corrections lift the smallness of Hj-iq^) to a certain extent, which makes it 
an observable that is particular sensitive to QED radiation. 

To supplement our calculation we carry out a dedicated Monte Carlo study, whose 
main purpose is three-fold. Eirst, we investigate how the electromagnetic logarithms are 
treated correctly in the presence of angular and energy cuts. We find that our analytical 
predictions can be directly applied, with the exception of the electron channel at BaBar, 
where our numbers have to be modified according to eqs. (7.1) and (7.2). Second, the size 
of the QED corrections, in particular their large relative size in HT{q^), are confirmed by 
the Monte Carlo (cf. tables 2 and 4). Last but not least, it consitutes also a validation 
of PHOTOS, which is used by experiments to estimate QED effects in the calculation of 
efficiencies. 

We update the Standard Model predictions for all angular observables integrated over 
two bins in the low-cp' region. The branching ratio and the observable 7Z(so) are also 
evaluated in the high-g^ region. Moreover, we provide our prediction for the zero crossing of 
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the forward-backward asymmetry (or, equivalently, 1-La)- The parametric and perturbative 
uncertainties are in general in the 5 — 15% range, exceptions are 1-La\^, 6] and the high-g^ 
branching ratio, where the relative errors are much larger. In the former case the reason 
is the zero crossing of Ha which entails a cancellation between the central values of the 
two bins in the low-(;^ region. In the latter case we suffer from poorly known hadronic 
parameters in the 1/m^’^ power-corrections, a drawback that is circumvented in the ratio 
lZ{so), which normalizes the B -A- Xsi^i~ rate to the inclusive -A Xuii^ rate with the 
same cut in [60]. 

We also study the sensitivity of the B -A Xsi'^l~ decay to new physics in a model- 
independent way. We give all observables in terms of ratios 8 , 9,10 of high-scale Wilson 
coefficients, which we assume to be altered by the new interactions. We also study cor¬ 
relations between different observables, bins and channels in the Rg — Riq plane, and 
extrapolate to the final Belle II data set of 50 ab“^. We find that Rt and Ha give the 
tightest constraints. On the other hand, if deviations from the Standard Model are seen, 
all observables become crucial to pin down the structure of new physics. 

In view of the recent measurement by LHCb [20] which reports a value for Rx = 
—7- K~^^~)/BR{B~^ -A K^e'^e~) in the low-q^ region that is significantly dif¬ 
ferent from unity, one might wonder whether this sign of lepton non-universality could 
be traced back to logarithmically enhanced QED corrections. LHCb uses the PHOTOS 
Monte Carlo to eliminate the impact of collinear photon emissions from the final state 
electrons. Therefore, the corrections calculated in this paper do not seem to apply to the 
ratio Rk- Given that the agreement between PHOTOS and our analytical calculations is 
not perfect (see e.g. tables 2 and 4), it would be advisable to correct for photon radiation 
using data-driven methods that do not rely on PHOTOS. 
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A QED and QCD functions 

A.l QED functions for the double differential rate 

Here we list the polynomials that appear in the functions of the log-enhanced 

QED corrections to the double differential rate in eq. (3.28). 
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Pl{s,z) 

P2{S,Z) 

P3{S,Z) 

P4{S,Z) 

P5{S,Z) 


P6{S,Z) 

P7{S,Z) 

P8{S,Z) 

P9{S,Z) 

Pio{s,z) 


Pu{s,z) 

P12{S,Z) 

P13{S,Z) 

P14{S,Z) 

P15{S,Z) 

P 16 {S,Z) 

P17{S,Z) 


25^ {z^ + 6z^ + l) +s {llz"^ - 8z'^ -3) + {z‘^ -if , 

45^ + l) + 3s^ - l) - 4s - l) - 9^2 - 7 , 

+s (6z® + 37z^ - 36^2 - 7) + + 24^^ + 3 , 

s^ (z® - 4z® + 2z^ - 28^2 - 3) - 35^ (z® - 4z® + - Az^ - l) 

+4s {z^ - 5z^ + 3z^ + l) -2 (5z^ + 24^^ + 3) , 

(13/ - 56/ + 210^"^ - 112/ - 55) 

+5^ (-15/ + 31/ - 127/ + 149/ + 154) 

+35^ (5z® - 9z® + 55/ - 312;^ - 84) 

+s (-13/ + 65/ - 285/ + 355^2 + 262) 

-13/ + 37/ - 299^2 - 109 , 
s (/ — l) — — 1 , 

s^ (43/ + 106^2 + 27) + 24s (2/ - z^ - l) + 3 (z^ - if , 
s^ (-z^° + 3/ + 32z^ + 364/ + 289/ + 17) 

+s (3z^° - 19z® + 106^*^ + 102/ - 173^2 - 19) 

+2(-/ + 7/-9/ + / + 2) , 

2s^ (17/ + 183/ + 143/ + 9) + s^ (77z® + 922/ - 92/ - 842^^ - 65) 
+s2 (/ - l)^ (46z® + 889^"^ + 1030^2 + 87) 

+s (/ - l)^ (256/ + 483/ + 51) + (/ - l)^ (74/ + ll) , 

-s® (I3z® - 66/ + 1288/ + 2706/ + 283) 

+s^ (-26z^° + 173/ - 2504/ - 2098/ + 7690/ + 989) 

+s^ (-13 z^ 2 + 122z^° - 1190/ + 830z® + 8809/ - 7288z^ - 1270) 

+s2 (/ - l)^ (15/ - 18/ + 397/ + 3716/ + 706) 

-s (/ - l)^ (15/ + 19/ - 403/ - 143) + (/ - l)^ (13/ - 22/ + l) , 
s^ (5z^ + 3) + / - 1 , 

s^ (z^ - 6/ - 9/ - 2) - s (z^ -lf + z^-1, 
s^ (/ + 22/ + 9) + s^ (/ + 11/ - 33/ - 11) 

-s (2/ + 17/ - 24z^ + 5) + (/ - l)^ (z^ + 7) , 
s^ (3/ + 12z^ + l) + s^ (4/ + 15/ - 18z^ - l) 

+s (z^ - l)^ (7/ - l) - (z'^ - l)^ , 
s (/ - l) + + 1 , 

s^ (5/ + 24z^ + 3) + s^ (6/ + 37/ - 36^^ - 7) 

+s (/ - 1)^ (19/ + 5) + (/ - 1)^ , 

s^ (z^ - 3/ + 39z^ + 27) - 2s (/ - 3/ - 5z^ + 7) + / - 3/ + 7^^ - 5 , 
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Pisis, z) = {3z^ - 18z^ - 49) - - lOz^ + 9) + 3s (z® - - llz^ + 16) 

-2Vs (z^ - 5) (z^ - 1)^ + (z^ - 1)^ (3z^ - 7) , 

Pi 9 (s, z) = (37^"^ + 86z^ + 21) + 16s {2z^ - - l) + - l)^ , 

P2o(s, z) = s^ [z^ - 4z^ + 154z^ + 340^2 + 85) - 2s (z® - 4z^ - + 28^^ + 33) 

+ (z^ - if (z^ - 2z^ + 5) , 

P2i(s, z) = s^ (3z^ - 37z^ - 359^2 - 183) - 2s®/2 (/ - 35^^^ - 5^^ + 39) 

+s2 (6z® - 77^® - 613^"^ + 305^2 + 379) - 4s^/2 (z^ - l)^ (z^ - 27z^ - 34) 

+s {z^ - if {3z^ - 31z^ - 323^2 - 229) - 2^/s (z^ - l)^ {z^ - Uz'^ - 29) 

+ {z^ - l)^ (3z^ - lOz^ - 33) , 

P 22 is, z) = 2s* (I5z® + 153z^ + 113^2 + 7) + s^ (69z^ + 754z® - 132/ - 642^^ _ 49) 

+s2 (/ - 1)^ (42z® + 717/ + 742^2 + 63) + 5s {z^ - l)^ (40/ + 63/ + 7) 

+ (/-l)^(38/ + 7) . (A.l) 


A.2 Functions for the QCD corrections to the T-Lj 

The one-loop QCD functions [49, 76] can be computed analytically, 


u 


( 1 ) 
77,T 


(/ 



f f^b'\ ' 

^/l -t- 1^ 

|2 ^ 

(^//2 lOs + 13 /j - 8^ 

1 Li2(l — s) 

\mbj 



6 (s - 1)2 



+ 

+ 


2 '/§ (/ — 6s — 3 ) Li 2 ^1 — '/s'j / ^ 3 //^ + 22s + 23/1 -|- 16^ ^/S — 1^ 


3(s- 1)2 
5/ - 5452 + 57s - 8 
18(s - 1)2 


36(s - 1)2 


- log(l - s) + 


s(5s + 1) log(s) 
3(s - 1)2 


+ -log(l - s) log(/ , 


'’79!r( 


4 

3 


/ /ifo \ 2v/(s + 3)Li2 1 

( 1 -/ 5 ) 

1 7r2 1 

(^16s 29/1 -h 19 ) (^/l - 1 ) 

\mb) 3(s — 1)2 


36(s - 1)2 



+ 


+ 


s-hl 8s- 


- s 


6(s — 1)2 


/ - 6s + 5 
6(s - 1)2 

(5s-hl)log(l-s) /3s-M)log(s) ,2 
-Si-+ 6(i-ip + 3 - ’> • 


CJ 


( 1 ) 
99,T 


(/ 


^/j + 1^ ^8s^/^ — 15s + 4/5 — 5^ Li2(l — s) 



2 (/ - 12s - 5) Li2 (1 - 
3(s — 1)2/| 

— 7s — 5) log(s) 

3(/-l)2 



- 43 - 


















s2 + 18s-19 (2 s + l)log(l-s) 2 


6(s — 1)2 


+ -log(l-s)log(s), 


(1) - 4 . / (2g^ -9^-3) Li,(l - g) 

710,a( ) 3 3(s - 1)2 3(,s - 1)2 

(3s 2 — 16s + 13) log ^1 — (4 s 2 — 13s — l) log ^1 — log(s) 

3(s - 1)2 3(s - 1)2 

(2s2 — 9s — 3) log(l — s) log(s) ^ (s^ — 23 s 2 + 23s — l) log(l — s) 

3(s - 1)2 6(s - l)2s 

fs - 20V§ + 5 ) f\/5 - l) ,^2 


6(s - 1)2 


(1) _ 2 (s2 - 3s - 1) Li 2 (l - s) 4(5 - 2s)sLi2 (l - Vtj [4:V§ - 3^ (^V§ - 1 

^910,a(s) - 3(s-1)2 3(s-1)2 3(s-1)2 

2 (2s2 — 7s + 5 ) log ^1 — 2 (s^ — 3s — 1 ) log(l — s) log(s) 

3(s- 1)2 3(s- 1)2 

(2s3 - lls2 + 10s - 1) log(l - s) ^ 2s(2s-5)log(l-\/i)log(s) ^2 

^ 3(s-l)2s ^ 3(s-l)2 y ’ 




(^+ 1 ) (4s^/^-7s + 2^/i-3jLi2(l-s) 9^2 _ 33 ^ 29 

“ "3 V^y ^ 3 (s-1)2\/5 6(s-1)2 


4 (s^ — 6s — 3 ) Li 2 ^1 — 
3(s - 1)2\/I 


+ 13s + 2^/§ + 9j — 1 

i8(s - iyv § 


(s^ - 3s + 2 ) log(l - s) 2 (s^ - 3s - 3 ) log(s) ^2 

-3(i3i)y-+ —3(i3ip— + 3 ’ 


4 f ^^b\ 4\/5(s + 3 )Li 2 (^1 - (^\/5 + l^ (^4s - 9\/| +3^ Li 2 (l - s) 

^79.l(« 1 = -3 log (^ —) + ^(-31)2 + ^(—1)2 

7s2-2s- 5 y (ss + 19\/5 + 5 ) - 1 ) (2s + l)log(l-s) 

6 (s- 1)2 18(s- 1)2 y 

(s —7)slog(s) 2 „ 

+ Hs-iy- +3Mi-»)i°gO). 
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^ — 8s + 3\/S — 4^ Li 2 (l — s) 4\/3 (s^ — 12s — 5) Li 2 — 


=- 


+ 


TT 


3(s-l)2 3(s-l)2 

(^35^/2 + 20s + \/5 + 8^ (^\/i - 4^3 _ ^ 1^2 _^^2s + 5 


+ 


18(s - 1)2 

8s(2s + l)log(s) 2 

3(s-l)2 + 3 Ml-^)log(.). 


+ 


6(s - 1)2 


- log(l - s) 


(A.2) 


The two-loop QCD functions [87, 88] are obtained from least-squares fits and are also valid 
for all q"^. The necessary data was kindly provided by the authors of [87, 88]. 

LUgfri^) = log (^) “*■ 54.919(1 - s)^ - 136.374(1 - s)^ 

-F 119.344(1 - s)2 - 15.6175(1 - s) - 31.1706, 


4wJ^) = log (^) + 74.3717(1 - s)^ - 183.885(1 - s)^ 

\Tnb / 

+ 158.739(1 - s)2 - 29.0124(1 - s) - 30.8056, 


= 4"^ log (—) - 5.95974(1 - s)^ + 11.7493(1 - s)^ 

-F 12.2293(1 - s) - 38.6457. (A.3) 

They are given for rih = 2 and n/ = 3. = 23/3 denotes the one-loop QCD /3-function 

for five active flavours. 

A.3 Functions for the QED corrections to the 77/ 

The following functions are again obtained by least-squares fits. They are valid in the 
low-g^ region (1 GeV^ < < 6 GeV^) only. 


, ,(em) 
^77 T 


is) 


, ,(em) 


(s) 


(em) 

‘^gg.T 




, ,(em) 

^gg,L 





1.54986 - 1703.72 s^ 1653.38 s^ - 683.608 s^ -h 179.279 s2 - 35.5047s 

8(1 - s )2 

9.73761 -g 647.747 s^ - 642.637 s^ 276.839 s2 - 68.3562 s - 4^ 

4(1-s)2 ’ 

2.2596 157.984 s^ - 141.281 s^ -p 52.8914 s^ - 13.5377 s -P 

2s(l - s)2 ’ 

-0.768521 - 80.8068 s^ -P 70.0821 s^ - 21.2787 s^ -p 2.9335 s - 

■ (l“1^ ^ 
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^79,TV®/ ~ 


In M 

\mj 


ml\ 19.063 + 2158.03.5^ - 2062.92s3 + 830.53s2 - 186.12s + 
mj ) 8(1 — s)2 


^79,T 1®/ 


In ^ 

\mj 


ml\ -6.03641 - 896.643s^ + 807.349s3 - 278.559s2 + 47.6636s - 
mj ) 4(1 — s)2 


(em)... _ fml\ \ 21.5291 + 3044.94s^ - 2563.05s3 + 874.074s2 - 175.874s + 0^4^ 

^27,T(®)-ln(^^j [ 8(1-s)2 

2.49475 + 598.376s^ - 456.8315^ + 117.683s2 - 9.90525s - Mli^" 

-I- ?- - - 

8(1-s)2 


+ ^4S(®1 ln( 


5GeV/ ’ 


(em). . [-8.01684 - 1121.13s^ + 882.711s^ - 280.866s2 + 54.1943s - 0128988 

-27.T(®1=ln^^J [ -4(135^- 

. -2.14058 - 588.771s^ + 483.9975^ - 124.579s2 + 12.3282s + Mi4^‘ 

-L 7 - - - 

4(1-s)2 


+ ^ (®1 In ( 


5GeV; ’ 


(emV . fmj\ [ 4.54727 + 330.182s^ - 258.194s^ + 79.8713s^ - 19.6855s + 0-037i348 

49.t^(®1 = In j [-^ 

. 73.9149s^ - 61.1338 s 3 + 14.6517s2 - 0.102331s + 0.7100371 


+ ^n( 


2s(l - s)2 


5GeV/ ’ 


(emV . fmj\ [-2.27221 - 298.369s^ + 224.662s^ - 65.1375s^ + 11.5686s - 0 0233098 

49,T^(®1 = In (y| j [-^ 

. -0.666157 - 120.303s^ + 109.3155^ - 28.2734s2 + 2.44527s + ^02797811 


(1 - s)2 


+ Y (®) In ( 


5GeV/ ’ 


4^4®^= In 3 


mj\ [2.84257 + 269.974s^ - 194.443s3 + 48.4535s2 - 8.24929s + 


2s(l - s)2 
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5GeV 


4(4.54727 + 330.182.s^ - 258.194s3 + 79.8713s2 - 19.6855s + 


0.0371348' 


9s(l - s)2 


64 (em).- , 2 / l^b 

+ 8l‘^99,T(^)ln ^ — 


5GeV; ’ 


(em). ^ fm?\ -1.71832 - 234.11s^ + 162.126s-^- 37.236ir +6.29949s- 


0.00810233 


5GeV 


8(224.662s3 - 2.27221 - 298.369s^ - 65.1375s9 + 11.5686s - ^ 


9(1 - s)2 


+ ^4“l’(«)1"^( ” 


5GeV7 ’ 


(em) /-.N _ 1 f 7-16\/I + 9s l + 3s M + \/5 \ s In s 


-riorA(«i= In 3 


4 (1 - s) 


+ ln(l - Vs) + 


(1-s) ’ 




5 —16\/5 + lls 1 —5 s M+\/5\ (1 —3 s) Ins 


4 (1-s) 


(1-s) 


(em) , fml\ r-351.322s^ + 378.173s3 - 160.158s2 + 24.2096s + 0.305176 


^210.AV®/ ~ ^n 2 


24s(l - s) 


7.98625 + 238.507 (s-a)- 766.869 (s-a)2.. , 2 ./. 

^^ -«) 


+ ln( 


5GeV 


with a = . 

The respective high-g^ functions for the branching ratio that are obtained by a least- 
squares fit (for fixed values of mj, and me) read 


r S 4 (s)+iS|(s) 1 16 (em). ^ , / hb 

^^29 (^) - ^n 2 , osi + n ^1010(«) ml 


mj) [8(1 -s)2(1 + 2s)J 9 


5GeVy ’ 


^5(g) ^41^; , / Pb 

22 Vm^y [8(1 -s)2(1 + 2s) 9(1-s)2(1 + 2s) V5GeV 


64 (em),.. , 2/ l^b 

+ ^^ioio(") 


5GeVy ’ 


33 (s) = In 


ml\ r S6(s) -HiS^(s) 1 8 

(em) ^ 7 hb 

3/ [ 96(1-s)2 +9‘^79 UgcV 


The functions Sj are polynomials in 5 = 1 — s and are valid for s > 0.65. 

S 4 (s) = -153.673 52 + 498.823 - 1146.74 + 1138.81 <5^ 
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(A.8) 


Si(s) = -255.712 < 5 ^ + 1139.10 - 2414.21 <5^ + 2379.91 <5^ , 

S 5 (.s) = -220.1015^ + 875.703 5^- 1920.56(5^ + 1822.07(5®, 

= -310.1135^ + 834.253 5^ - 2181.945^ + 2133.785® , 

S^(s) = -518.180 52 + 2047.18 5^-4470.04 5^ + 4827.74 5® . 

B New Physics formulas 

nTil, 3.5]ee = [0.0162226X(i?7i?^) + 0.00186782X(ii7i?3) + 0.00985919X(i?8^9) 

- 0.000201564X(7?8-Rto) + 0.0465868X(i?7) - 0.00822885X(ii8) 

- 0.0187815X(i?9) + 0.000379966X(i?io) + 0.3931567^(i^7) 

+ 0.0400072 n{Rs) + 0.0531851 TZiRg) - 0.0385002 7^(i^lo) 

+ 0.04584277^(i^7i^8) - 0.3699647^(i^7i^;) + 0.005706077^(i^7i^lo) 

- 0.03694987^(i^8-R9) + 0.0006164227^(i^8^^o) “ 0.00978058 7^(i^9i^^o) 

+ 0.204994 |i?7|2 + 0.00230146 + 0.244813 |i?9|2 

+ 1.74294 |7?io| 2+ 0.632156] x 10■^ (B.l) 

'HT[3.5,6]ee =[0.00519889X(ii7i7^) +0.00141211 X(i?7i?^) + 0.00745377X(i?8^9) 

- 0.000152386X(i?8-Rto) + 0.0151043X(ii7) + 0.00358335X(ii8) 

- 0.0100672X(i?9) + 0.000148662X(i7io) - 0.1385167^(i^7) 

- 0.0131665 7^(7^8) + 0.375959 7^(7^9) - 0.074623 7^(i^lo) 

+ 0.01435687^(i^7i^8) - 0.254325 7^(i^7i?;) + 0.004311397^(i^7i?lo) 

- 0.02609437^(i^8-R9) + 0.0004676877^(i^8^^o) “ 0.01572597^(i^9i^^o) 

+ 0.0631028 |i?7|2 + 0.000727107 |i?8p + 0.273706 |i?9|2 
+ 1.96638 |i?io|2+ 0 . 257773 ] X 10■^ (B.2) 

6 ]ee = 0.0214215X(i?7i?^) + 0.00327993X(ii7ii3) + 0.017313X(i?87?9) 

- 0.000353949X(7?8-Rto) + 0.0616911 X(ii7) - 0.0046455X(i78) 

- 0.0288487X(i?9) + 0.000528628X(i?io) + 0.254647^(i^7) 

+ 0.02684077^(i^8) + 0.429144 7^(i^9) - 0.113123 7^(i^lo) 

+ 0.06019947^(i^7i^^) - 0.6242897^(i^7i^;) + 0.01001757^(i^7i^^o) 

- 0.0630441 niRsRl) + 0.00108411 n{RsRlo) “ 0.0255065 
+ 0.268097 |i?7|2 + 0.00302857 |i? 8 p + 0.518519 \Rg\‘^ 

+ 3.70932 |7?ioP+ 0 . 889929 ] x 10■^ (B.3) 

'Ht[ 1,3.5]^/, = 0.0162226X(i?7i?^) + 0.00186782X(i?7i?3) + 0.00985919X(i?87?;) 

- 0.000201564X(i?8-Rto) + 0.0478295X(ii7) - 0.00813434X(i?8) 
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^A[l,3.5]ee 


^A[3.5,6]ee 


- 0.0247652X(i?9) + 0.000379966X(i?io) + 0.45956377(i?7) 

+ 0.0451794 77(7?8) - 0.155638 77(7?9) - 0.0385002 77(i?io) 

+ 0.046052177(i?7i?^) - 0.33743177(i?7i?;) + 0.0057060777(i?7i?io) 

- 0.034475777(i?8-R9) + 0.00061642277(i?8^to) “ 0.00978058 77(i?9-Rto) 

+ 0.206371 |i?7|2 + 0.00230943 |i?8p + 0.179467 |i?9p 


+ 1.288811 + 0.436438 


X 10 “^ 


(B.4) 


0.00519889 X(ii7ii8) + 0.00141211 X(i?77?^) + 0.00745377 X(i?87?9) 


- 0.000152386X(i?8-Rto) + 0.0165184X(ii7) + 0.00369089X(i?8) 

- 0.0169196X(i? 9 ) + 0.000148662X(i?io) - 0.11237677(7?7) 

- 0.0111424 77(i?8) + 0.24902777(i?9) - 0.074623 77(i?io) 

+ 0.014654777(i?7i?8) - 0.24467177(i?7i?;) + 0.0043113977(i?7i?io) 

- 0.0253601 77(7?8-R9) + 0.00046768777(i?8^to) “ 0.015725977(i?9i?io) 

+ 0.0650616 |i?7p + 0.000738436 |i?8p + 0.239011 


+ 1.72527 |7?ioP + 0.123204 


X 10 ■^ 


(B.5) 


0.0214215X(i?7i?^) + 0.00327993X(ii7i?3) + 0.017313X(i?8^9) 


- 0.000353949X(i?8-Rto) + 0.0643479X(ii7) - 0.00444346X(i?8) 

- 0.0416848X(i? 9 ) + 0.000528628X(i?io) + 0.34718677(i?7) 

+ 0.034037 77(i?8) + 0.0933889 77(7?9) - 0.113123 77(i?io) 

+ 0.0607068 niRrRl) - 0.582101 n{RjRl) + 0.0100175 n{R^RlQ) 

- 0.0598358 TZiRsRg) + 0.00108411 n{RsRlo) “ 0.0255065 TZiRgRlo) 
+ 0.271433 |i?7p + 0.00304786 |i?8p + 0.418478 |i?9p 


+ 3.01408 |i?ioP + 0.559642 


X 10 ■^ 


(B.6) 


- 0.0000761415X(ii8i?9) + 0.0259112X(ii8i?io) + 0.0031943X(ii9i?^o) 


- 0.000083788X(7?8) + 0.00025712X(7?9) - 0.112552X(i?io) 

+ 0.023027777(i?7) + 0.0018154377(i?8) - 0.0133235 77(ii9) 

- 0.82662677(iiio) + 0.00214715 77 (i? 7 ^ 9 ) - 0.84915477(i?7i?^o) 

+ 0.00022240177(i?8^9) - 0.084738977(i?8^to) + 0.72293477(i?9i?^o) 


-0.00174093 Iii9p 


0.0120987 |i?ioP + 0.0121072 


X 10 ■^ 


(B.7) 


- 0.000057133X(i?8^3) + 0.0194427X(i?8i?l:o) + 0.00509883X(ii9i?^o) 
- 0.000062727X(7?8) + 0.000151953X(i?9) - 0.0912157X(7?io) 

+ 0.0172872 77(i?7) +0.0013674477(i?8) - 0.0259495 77(ii9) 
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+ 0.3562937^(i^lo) + 0.001603797^(i^7^9) “ 0.6051037^(i^7-R^o) 

+ 0.000169807- 0.0623319+ 1-08406 

- 0.0027675 |i?9p-0.0192329 liiiol^- 0.0115297] x 10■^ (B.8) 

UaII, 6]ee = \ - 0.000133274X(i?8-R9) + 0.0453539X(i?8-Rl:o) + 0.00829314X(ii9i?^o) 

- 0.000146515X(7?8) + 0.000409073X(i?9) - 0.203767X(iiio) 

+ 0.04031497^(i^7) + 0.003182877^(i^8) - 0.03927317^(i^9) 

- 0.4703337^(i^lo) + 0.003750947^(i^7i^;) - 1.454267^(i^7i^lo) 

+ 0.0003922097^(i^8-R9) - 0.1470717^(i^8-R^o) + 1-806997^(i^9i^^o) 

- 0.00450843 Ii?9p-0.0313316 I+ 0.000577448] x 10■^ (B.9) 

HaII, 3.5]^^ = f - 0.0000761415X(ii8i79) + 0.0259112X(i?8^1:o) + 0.0031943X(ii9i7^o) 

- 0.000083788X(7?8) + 0.00025712X(7?9) - 0.112552X(i?io) 

+ 0.02302777^(i^7) + 0.001815437^(i^8) - 0.0133235 7^(i^9) 

- 0.8756077^(i^lo) + 0.00214715 7^(i^7i^3) - 0.8453277^(i^7i^^o) 

+ 0.0002224017^(i^8-R9) - 0.08444787^(i^8-R^o) + 0.694542 7^(i^9i?^o) 

- 0.00174093 |ii9p-0.0120987 |i?ioP+ 0.0131242 x 10■^ (B.IO) 

'Ha[ 3.5, 6]^^ = [ - 0.000057133X(i?8-R3) + 0.0194427X(i?8-Rl:o) + 0.00509883X(ii9i?^o) 

- 0.000062727X(7?8) + 0.000151953X(i?9) - 0.091289X(iiio) 

+ 0.01728727^(i^7) +0.001367447^(i^8) - 0.0259495 7^(i^9) 

+ 0.3180087^(i^lo) + 0.001603797^(i^7i^9) - 0.619516 

+ 0.0001698077^(i^8-R9) - 0.0634287^(i^8-R^o) + 1-077867^(i^9i^^o) 

- 0.0027675 |i?9p-0.0192329 |i?io|^- 0.0113078 x 10■^ (B.ll) 

'Ha[ 1, 6]^^ = [ - 0.000133274X(i78-R9) + 0.0453539X(i78-Rl:o) + 0.00829314X(ii9i7^o) 

- 0.000146515X(7?8) + 0.000409073X(i?9) - 0.203841 X(iiio) 

+ 0.04031497^(i^7) + 0.003182877^(i^8) - 0.03927317^(i^9) 

- 0.5575997^(i^lo) + 0.003750947^(i^7i^;) - 1.464847^(i^7i^^o) 

+ 0.0003922097^(i^8-R9) - 0.1478767^(i^8-R^o) + 1-77241 

- 0.00450843 liigp-0.0313316 |i?ioP+ 0.00181642 x 10■^ (B.12) 

'H3[l,3.5]ee = 0.0264036X(i?io) + 3.071567^(i^lo) - 1.740437^(i^7i^lo) 

- 0.132357 7^(i^8i^^o)+ 2.94364 7^(i^9i^^o)-0.105444] x 10■^ (B.13) 

-^3(3.5, 6]ee = 0.132813X(i?io) + 3.519047^(i?lo) - 0.9133537^(i^7i^l;o) 
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^3[l,6]ee 

^3[1,3.5]^^ 

^3 [3.5,6]^^ 

^3 [1,6]^;. 

^4[l,3.5]ee 

^4[3.5,6]ee 

^4[l,6]ee 

^4[1,3.5]^^ 


- 0.0694587 + 2.4359 77(i?9i?io) - 0.0872558 


X 10 


-9 


(B.14) 


p. 159216X(i?io) + 6.590677(i?io) - 2.6537977(i?7i?io) 
- 0.201815 77(i?8^1:o) + 5.37954 77(i?9i?^o) “ 0.192699 


X 10 


-9 


(B.15) 


p.010976X(i?io) + 1.2794677(i?io) - 0.723502 77(i?7i?^o) 
- 0.0550209 77(i?8.Rto) + 1.22368 77(ii9i?io) “ 0.0438331 


X 10■^ 

(B.16) 


0.0552105X(i?io) + 1.46503 77(i?io) - 0.379682 77(i?7i?^o) 
-0.028874177(i?8i?w) + 1-0126177(i?9i?l;o) - 0.03627241 x 10“® 


(B.17) 


p.0661865X(i?io) + 2.7444977(i?io) - 1.1031877(i?7i?io) 
- 0.083895 77(i?8^1:o) + 2.23628 77(i?9i?io) - 0.0801055] 


X 10 


-9 


(B.18) 


- 0.0412679X(i?7) - 0.00313835X(i?8) + 0.200198X(i?9) 


- 0.43003477(ii7) - 0.03405877(ii8) + 1.4651677(i?9) 

0.0135748 77(i?7i?^) - 0.361104 77(i?7i?;) - 0.0274613 77(i?8.R9) 

0.482688 \Rg\‘^ + 0.0892516 |i?7p + 0.00051617 |i?8p 


3.35446 17?ioP + 1.6742 


X 10 


-9 


(B.19) 


- 0.0257056X(i?7) - 0.00195486X(i?8) + 0.127314X(i?9) 


- 0.17595 77(i?7) - 0.0138586 77(ii8) + 0.528054 77(ii9) 

0.00348411 TZiRjRl) - 0.127392 TZiR^Rl) - 0.00968792 TZiRsRl) 
+ 0.179914 |i?9|2 + 0.0229073 |i?7p + 0.00013248 |i?8p 


1.25032 |i?ioP + 0.529364 


X 10 


-9 


(B.20) 


- 0.0669735 X(7?7) - 0.0050932X(i?8) + 0.327512 X(7?9) 


- 0.605984 77(i?7) - 0.0479166 77(7?8) + 1.99322 77(ii9) 

+ 0.0170589 77(i?7i?8) - 0.48849677(i?7i?;) - 0.037149277(i?8.R9) 

+ 0.662601 \Rg\‘^ + 0.112159 {RjI'^ + 0.00064865 [iisp 


+ 4.60478 |7?ioP + 2.20357 


X 10■^ 


(B.21) 


^ - 0.0171551 X(ii7) - 0.00130462 X(i?8) + 0.0832226X(ii9) 

- 0.179086 77(ii7) - 0.0141823 77(i?8) + 0.609926 77(i?9) 

+ 0.0056430877(i?77?8) - 0.150112 77(i?7i23) - 0.011415777(ii87?9) 
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(B.22) 


+ 0.200654 \Rgf + 0.0371021 |i?7p + 0.000214573 
+ 1.39446 |7?ioP+ 0.697498] x 10■^ 

'H4[3.5, = - 0.0106858X(i?7) - 0.000812638X(i?8) + 0.0529245X(7?9) 

- 0.073255777(i?7) - 0.0057696477(i?8) + 0.219832 77(729) 

+ 0.0014483577(i?7728) - 0.052957177(72772^) - 0.0040272977(72872^) 

+ 0.0747905 | 729p + 0.00952261 |727p + 0.0000550722 | 728p 

+ 0.51976 |72ioP+ 0.22061 x 10■^ (B.23) 

= [- 0.027841 X(727) - 0.00211725X(728) + 0.136147X(729) 

- 0.25234177(727) - 0.0199519 77(728) + 0.829758 77(729) 

+ 0.0070914377(727728) - 0.20306977(727729) - 0.01544377(728729) 

+ 0.275445 |729|2 + 0.0466247 |727p + 0.000269645 |728|^ 

+ 1.91422 |72ioP+ 0.918108] x 10■^ (B.24) 

77L[l,3.5]ee =[o.000741931X(72772^) + 0.000952641X(727723) +0.0050284X(728723) 

- 0.000102803X(72872to) + 0.00124959X(727) + 0.00594309X(728) 

+ 0.00735758X(729) - 0.00113202X(72io) - 0.19486677(727) 

- 0.0251935 77(728) + 1.4250177(729) - 0.25154 77(72io) 

+ 0.0021375177(727728) - 0.13628377(727729) + 0.00300802 77(72772io) 

- 0.0187453 77(72872^) + 0.00040242177(72872io) - 0.046284177(72972io) 

+ 0.00589466 |727p + 0.000128527 | 728p + 0.575967 | 729p 

+ 4.20578 |72ioP+ 0.806915] x 10■^ (B.25) 

77L[3.5,6]ee =f0.000562052X(72772^) + 0.000724099 X(727723) +0.00382208 X(72872;) 

- 0.0000781401 X(72872io) + 0.00223749 X(727) + 0.0047901 X(728) 

- 0.00211229X(729) - 0.000740423X(72io) - 0.16111777(727) 

- 0.019209477(728) + 1.14892 77(729) - 0.193345 77(72io) 

+ 0.001762477(727728) - 0.10750177(72772^) + 0.0022863677(72772io) 

- 0.0136079 77(72872^) + 0.00028642377(72872^9) - 0.035510977(72972);o) 

+ 0.00631092 |727p + 0.0000975709 | 728p + 0.439598 | 729p 

+ 3.20293 |72ioP+ 0 . 701014 ] x 10■^ (B.26) 

77^[1, 6]ee = 0.00130398X(727728) + 0.00167674X(727729) + 0.00885049X(728729) 

- 0.000180943X(72872to) + 0.00348707X(727) + 0.0107332X(728) 
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+ 0.00524529X(i?9) - 0.00187244X(i?io) - 0.3559827^(i^7) 

- 0.0444029 7^(4^8) + 2.57393 7^(i^9) - 0.444885 7^(i^lo) 

+ 0.003899917^(i?7i^8) - 0.2437847^(i^7i^3) + 0.005294387^(i^7i^lo) 

- 0.0323532 7^(i^8-R9) + 0.0006888437^(i^87^^o) “ 0.081795 7^(i^9i?l;o) 

+ 0.0122056 |i?7p + 0.000226098 |i?8p + 1.01556 |i?9p 
+ 7.40871 |i?ioP + 1.50793 x 10■^ (B.27) 

'Hl[ 1,3.5]^^ = 0.000741931 X(7?7i?8) + 0.000952641 X(i?7i?3) +0.0050284X(ii8i?9) 

- 0.000102803X(i?8.Rto) + 0.000345511 X(i?7) + 0.00587433X(i?8) 

+ 0.0117155X(i?9) - 0.00113202X(i?io) - 0.2172457^(i^7) 

- 0.0269584 7^(i^8) + 1.53068 7^(i^9) - 0.25154 7^(i^lo) 

+ 0.002605737^(i?7i^8) - 0.1530577^(i^7i^9) + 0.003008027^(i^7i^^o) 

- 0.0200209 7^(7^8.R9) + 0.0004024217^(i^8.R^o) “ 0.04628417^(i^9i^^o) 

+ 0.00897313 |i?7p + 0.000146331 \Rs\^ + 0.609248 |i?9p 
+ 4.43707 |7?ioP+ 0.914888 x 10■^ (B.28) 

= 0.0005620521{Rj Rl) + 0.0007240991{RjRl) + 0.003822081(RsRg) 

- 0.0000781401 X(i?87?to) + 0.00132426X(i77) + 0.00472065X(i?8) 

+ 0.00235817X(i?9) - 0.000740423X(i?io) - 0.1773927^(i77) 

- 0.0204867 7^(i^8) + 1.239117^(i^9) - 0.193345 7^(i^lo) 

+ 0.001940947^(i?7i^^) - 0.1180587^(i^7i^9) + 0.002286367^(i^7i^lo) 

- 0.0144108 7^(i^8.R9) + 0.0002864237^(i^8.R^o) “ 0.03551097^(i^9i^lo) 

+ 0.00748476 |ii7p + 0.00010436 |i?8p + 0.466941 {Rgl"^ 

+ 3.39295 |i?ioP+ 0.791074 x 10■^ (B.29) 

77l[ 1,6]^^ = 0.00130398 X(ii7i?^) +0.00167674X(i?77?^) + 0.00885049 X(i?8.R9) 

- 0.000180943X(7?87?to) + 0.00166977X(7?7) + 0.010595X(i78) 

+ 0.0140737X(i?9) - 0.00187244X(i?io) - 0.3946387^(i^7) 

- 0.04744517^(i^8) + 2.76979n{Rg) - 0.444885 7^(i^lo) 

+ 0.004546677^(i?7i^8) - 0.271115 7^(i^7i^9) + 0.005294387^(i^7i^lo) 

- 0.03443177^(i^8.R9) + 0.0006888437^(i^8.R^o) “ 0.081795 7^(i^9i^lo) 

+ 0.0164579 |i?7p + 0.000250691 |i?8p + 1.07619 |i?9p 
+ 7.83003 |7?ioP + 1.70596 x 10■^ (B.30) 

B[l, 3.5]ee = [0.0169646X(i?7i?^) + 0.00282046X(ii7i?3) + 0.0148876X(i?8-R9) 
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- 0.000304367X(i?8-Rto) + 0.0347138- 0.00283044X(i?8) 

+ 0.000660238X(7?9) - 0.00100106X(7?io) + 0.189792 7^(i?7) 

+ 0.0139496 n{Rs) + 1.462717^(i^9) - 0.290285 7^(i^lo) 

+ 0.0507378 7^(i^7-R8) “ 0.5062517^(i^77^;) + 0.008714097^(i^7i?^o) 

- 0.0584716 7^(i^8-R9) + 0.001076437^(i^8-R^o) “ 0.05606477^(i^9i?lo) 

+ 0.210889 |i?7|2 + 0.0028916 |i78p + 0.813297 

+ 5.94874 |i?ioP + 1.46402] x 10■^ (B.31) 

.B[3.5,6]ee =[0.00576094X(i?7i7^) +0.00213621 X(i?77?^) + 0.0112758X(i?8^9) 

- 0.000230526X(7?8-Rto) + 0.0117001 X(i?7) + 0.00792519X(ii8) 

- 0.000973809X(i? 9 ) - 0.000822616X(i?io) - 0.3041977^(i^7) 

- 0.0338418 7^(i^8) + 1.5387^(i?9) - 0.2682057^(i?lo) 

+ 0.0166482 - 0.361825 7^(i^7i^;) + 0.00659775 7^(i^7i^lo) 

- 0.0407383 7^(i^8-R9) + 0.000775603“ 0.05123687^(i^9i^lo) 

+ 0.0694138 |i?7p + 0.000881518 |i?8p + 0.714084 {Rgl"^ 

+ 5.16931 |i?ioP +0.985134] x 10■^ (B.32) 

B[l, 6]ee = [0.0227255X(i?7i?^) + 0.00495667X(ii7i?3) + 0.0261634X(i?8-R9) 

- 0.000534893X(i?8-Rto) + 0.0464139X(ii7) + 0.00509475X(i?8) 

- 0.000313571 X(7?9) - 0.00182368X(i?io) - 0.114406 7^(i?7) 

- 0.0198921 n{Rs) + 3.000717^(i^9) - 0.55849 7^(i^lo) 

+ 0.0673867^(i^7i^8) - 0.868076+ 0.01531187^(i^7i^lo) 

- 0.0992099 7^(i^8-R9) + 0.001852037^(i^8-R^o) “ 0.1073017^(i^9i^^o) 

+ 0.280302 |i?7|2 + 0.00377311 jiisp + 1.52738 |i?9p 
+ 11.1181 |i?ioP + 2 . 44915 ] X 10■^ (B.33) 

B[l, 3.5]^^ = 0.0169646X(i?7i?^) + 0.00282046X(ii7i?3) + 0.0148876X(i?8-R9) 

- 0.000304367X(7?8-Rto) + 0.0350544X(i?7) - 0.00280454X(ii8) 

- 0.000975567X(i?9) - 0.00100106X(i?io) + 0.2338327^(i?7) 

+ 0.017358 niRs) + 1.35952 7^(i^9) - 0.290285 7^(i^lo) 

+ 0.0514155 7^(i^7i^8) - 0.4904897^(i^7i^;) + 0.008714097^(i^7i^lo) 

- 0.0572729 7^(7^8-R9) + 0.001076437^(7^8-R^o) “ 0.05606477^(i?9i^^o) 

+ 0.215344 IRjI'^ + 0.00291736 |i?8p + 0.78123 |i?9p 
+ 5.7259 |iiioP + 1.3762] x 10■^ (B.34) 

B[3.5,6]f,f, = 0.00576094X(i?7i7^) +0.00213621 X(i?77?^) + 0.0112758X(i?87?9) 
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- 0.000230526X(i?8-Rto) + 0.0122024X(ii7) + 0.00796339X(i?8) 

- 0.00336638X(ii 9 ) - 0.000822616X(7?io) - 0.2943377(7?7) 

- 0.0330948 77(i?8) + 1.50123 77(ii9) - 0.268205 77(i?io) 

+ 0.017124777(i?7-R8) “ 0.36272877(i?77?;) + 0.00659775 77(i?7i?w) 

- 0.04080777(i?8^9) + 0.00077560377(i?8-Rto) “ 0.051236877(i?9i?io) 

+ 0.0725464 |i?7p + 0.000899635 |i?8p + 0.706733 

+ 5.11822 I+ 0.940534] x 10■^ (B.35) 

.B[l,6]^^ =[0.0227255X(i?7i?^) + 0.00495667X(i?7ii3) + 0.0261634X(7?8-R9) 

- 0.000534893X(7?8-Rto) + 0.0472568X(i?7) + 0.00515885X(ii8) 

- 0.00434195X(i? 9 ) - 0.00182368X(i?io) - 0.0604983 77(i?7) 

- 0.0157368 77(i?8) + 2.86075 77(ii9) - 0.55849 77(i?io) 

+ 0.0685402 77(i?7i?^) - 0.85321777(i?7i?;) + 0.015311877(i?7i?io) 

- 0.0980877(i?8^9) + 0.0018520377(i?8^w) “ 0.10730177(i?9i?io) 

+ 0.287891 \R^\‘^ + 0.003817 |i?8|^ + 1-48796 {Rgl"^ 

+ 10.8441 |i?ioP +2.31673] x 10■^ (B.36) 

B[> UA]ee =f0.000264356X(i?7i?^) + 0.000401975X(ii7i?3) +0.00161219X(i?8-R9) 

- 0.0000328066X(i?8-Rto) “ 0.0158129X(i?7) + 0.000478008X(i?8) 

+ 0.125395X(7?9) - 0.00293188X(i?io) - 0.0723471 TZiRr) 

- 0.0082779377(i?8) + 0.51171577(i?9) - 0.0806142 77(i?io) 

+ 0.00070967877(i?7i?8) - 0.0516424+ 0.0011161477(i?7i?to) 

- 0.00651216 77(i?8-R9) + 0.000119004 77(i?8-Rl;o) “ 0.0168936 77(i?9ii^o) 

+ 0.00287361 |i?7p + 0.0000373632 |ii8p + 0.211548 {Rgl"^ 

+ 1.50748 |i?ioP+ 0.200589] x 10■^ (B.37) 

13[> 14.4]^^ = 0.000264356X(i?7i?^) + 0.000401975 X(ii7i?3) +0.00161219X(i?8-R9) 

- 0.0000328066X(i?8-Rto) “ 0.0175987X(7?7) + 0.000342205X(i?8) 

+ 0.134924X(i? 9 ) - 0.00293188X(i?io) - 0.087186377(ii7) 

- 0.0094385277(i?8) + 0.59439377(i?9) - 0.0806142 77(i?io) 

+ 0.00083552777(i?7i?8) - 0.060198477(i?7i?^) + 0.0011161477(i?77?to) 

- 0.00716282 77(i?87?^) + 0.000119004 77(i?8-Rw) “ 0.0168936 77(i?9ii^o) 

+ 0.00370104 |ii7|2 + 0.0000421485 |i?8p 

+ 0.234333 Ii ?9 p + 1.66583 Ii?ioP+ 0.292268] x 10■^ (B.38) 

77(14.4)ee = 0.000352294 X(7?7i?^) + 0.000544926 X(i?7i?3) + 0.00213997X(7?8-R9) 
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- 0.0000442492X(i?8^to) “ 0.0160419X(i?7) + 0.000523537X(i?8) 

+ 0.130938X(4?9) - 0.00323922X(7?io) - 0.066941177(i?7) 

- 0.00821459 n{R^) + 0.458105 n{R<^) - 0.095890177(i?io) 

+ 0.00080755877(i?77?8) “ 0.05486477(i?7i?;) + 0.0012343277(i?7ii^o) 

- 0.00734198 77(i?8^9) + 0.000139543 77(i?8i?w) “ 0.0189772 77(i?9i?^o) 

+ 0.00293717 |ii7p + 0.0000444449 jiisp + 0.228597 {Rgl"^ 


+ 1.6322 liiiol^ + 0.174573 


X 10■^ 


(B.39) 


77 ( 14 . 4 )^^ = 


0.000352294 X(i?7i?^) + 0.000544926 X(i?7i?3) + 0.00213997X(7?8-R;) 


- 0.0000442492X(i?87?to) “ 0.0181914X(i?7) + 0.000360068X(i?8) 

+ 0.142407X(i?9) - 0.00323922X(i?io) - 0.082515477(ii7) 

- 0.0094376277(i?8) + 0.5454477(i?9) - 0.095890177(i?io) 

+ 0.00095904477(i?7i?^) - 0.065163177(i?77?9) + 0.0012343277(i?7i?to) 

- 0.008125277(i?8-R9) + 0.00013954377(i?87?to) “ 0.0189772 77(ii9i?^o) 

+ 0.00393316 |ii7|2 + 0.000050205 |i?8p + 0.256024 \Rq\‘^ 


+ 1.822811 + 0.266662 


X 10“^ 


(B.40) 
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